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ABSTRACT 


A  theoretical  prediction  by  Dr.  Avihu  Ze'evi  of  the 
correlation  between  direct  path  and  folded  path  optical 
scintillation  was  developed  and  published  at  the  Naval  Post¬ 
graduate  School  in  March  1982.  Simultaneously,  a  method  for 
measuring  the  scintillation  of  two  laser  beams  propagating 
along  direct  and  folded  paths  was  devised  and  used  to  record 
preliminary  findings  supporting  the  theoretical  work.  This 
system  has  been  improved  upon  to  permit  a  complete  experi¬ 
mental  verification  of  the  theory  to  be  attempted  outside 
the  laboratory.  The  experimental  results  do  not  support 
Dr.  Ze'evi' s  theory;  however,  changing  turbulence  conditions 
while  measuring  the  scintillation  strength  may  have  been 
the  cause.  Further  attempts  at  verification  will  have  to 
control  this  factor. 
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I.  INTRODUCTION 


A.  BACKGROUND 

When  light  from  a  constant  source  propagates  through  the 
atmosphere,  the  intensity  distribution  at  a  point  on  the 
receiver  is  altered  by  inhomogeneities  in  the  atmospheric 
index  of  refraction  along  the  path.  In  simple,  quantitative 
models,  the  index  of  refraction  is  primarily  a  real  function 
of  temperature  and  pressure,  and  has  the  following  form: 

n(r)  =  nQ  +  n^r) 

where  n(r)  is  the  real  index  of  refraction 

at  a  point  r, 

nQ  *  <n(r,t)>  is  the  time  averaged 
(ensemble)  value  of  n(r)  at 
a  point  r,  usually  given  the 
value  of  1,  and 

nj(r)  is  the  perturbation  in  refractive 
index  at  a  point  r  such  that 
nj (r)  <<  nQ . 

If  we  assume  that  pressure  fluctuations  in  the  atmosphere 
are  small  (an  isobaric  atmosphere)  then  the  index  of 
refraction  at  a  point  becomes  a  function  of  temperature. 

Within  the  atmosphere  near  the  Earth's  surface,  point 
temperature  inhomogeneities  are  caused  by  convection  on  a 
large  scale  and  random  molecular  motion  on  the  microscale. 
Large  scale  vertical  temperature  gradients  in  this  boundary 
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layer  generate  large  scale  motion  of  air  masses,  imparting 
energy  to  the  velocity  field  and  creating  turbulent  eddies 
of  a  size  comparable  with  the  inversion  height  or  height 
above  ground.  In  the  Kolmogorov  model  of  turbulence,  these 
large  vortices  break  down  into  smaller  and  smaller  vortices 
until  the  flow  stabilizes  at  a  minimum  vortex  size  where 
the  local  Reynolds  number  within  the  vortex  is  less  than 
the  critical  value  for  air. [Ref.  1].  The  outer  scale  of 
the  turbulence  (diameter  of  the  large  scale  vortex)  is 
termed  LQ  and  the  inner  scale  (diameter  of  the  small  scale 
vortex)  is  termed  lQ.  Turbulent  vortices  throughout  the 
spectrum  of  sizes  from  LQ  to  result  in  variations  in 
the  index  of  refraction  at  any  given  point  in  the  atmosphere 
as  a  function  of  time. 

An  accepted  measure  of  the  turbulence  in  the  atmosphere 

is  the  refractive  index  structure  constant,  C  2  .  C  2  is 

n  n 

proportional  to  n^r)  and  may  be  derived  from  wire  probe 
temperature  difference  sensors.  Typical  values  for  Cn2  in 
the  lower  atmosphere  range  from  10~17  (m'2/3)  for  weak 
turbulence  to  10'12(ra'2/3)  for  strong  turbulence.  As  one 
would  expect  for  a  temperature  dependent  parameter  Cn2  , 
varies  diurnally  and  reaches  a  peak  during  the  midda^  period. 
[Ref.  2], 

Direct  path  optical  scintillation  may  be  defined  as 
fluctuations  in  the  irradiance  (common  usage  is  to  refer  to 
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the  irradiance  as  intensity)  at  the  receiver,  due  to  inhomo¬ 
geneities  in  the  index  of  refraction  over  the  propagation 
length  causing  random  dispersion  of  the  light  as  a  function 
of  space  and  time.  Because  Cn2  is  a  direct  measure  of  the 
level  of  atmospheric  turbulence,  it  also  varies  directly 
with  the  amount  of  scintillation  observed,  and  appears  in 
equations  describing  the  strength  of  the  scintillation. 

An  accepted  measure  of  the  scintillation  strength  is 
the  variance  of  the  normalized  intensity,  as  follows: 

CTj2*  <(I/IQ)2  >  -  <  I/IQ> 2 

T~ 

°«  [<I2>  -  <I>21/I02 

■  °i  '  K 

where  I  =  I(r,t)  is  the  intensity  on  the 

receiver  at  an  instant  in  time, 

I  =  <I(r,t)>  is  the  ensemble  average, 

0  and 

CJj2  =  <I2>  -  <I>2. 

Experimentally  it  has  been  found  that  if  we  restrict  ourselves 
to  regions  of  weak  turbulence,  then  the  log  intensity  of 
the  electromagnetic  field  is  also  a  random  variable  with 
normal  distribution  [Ref.  3],  It  is  therefore  convenient 
to  describe  the  scintillation  strength  as  the  log  intensity 
variance,  as  follows: 
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then 

and 


Thus , 


£  =  log  (I/I0) 

a  2  *  <l2>  -  <l>2 
l 

0£2  =  <  [  log  (y  )2]>  -  <l0g  (y  )>2 
o  o 

=  <(log  I  -  log  IQ) 2 >  -  <log  I  -  log  IQ>2 
*  <(log  I)2>  -  2  log  IQ  <log  I>  + 

(log  I0)2  -  < log  I>2  -  (log  IQ)2 
+  2  log  IQ  < log  I> 

=  < (log  I)2>  -  <log  I>2 


a 


2 

logl. 


In  this  way  the  variance  of  the  received  signal  intensity  is 
seen  to  have  no  dependence  upon  the  average  signal  intensity, 
and  the  log  normal  distribution  may  be  completely  described 
by  a  single  parameter,  a  2  . 

Assuming  conditions  of  homogeneous  turbulence  with 
5-0  <<  (*  L)  A  ,  the  log  amplitude  variance  was  shown  by 

Tatarski  [Refs.  3,  4]  to  obey  the  following  relations: 
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Plane  Waves 

ax2  -  0.31  K?/6 

C  2 

Ln  L 

Spherical  Waves 

ax2  -  .125  K  7/6 

where 

K  *  2tt/X , 

L  =  path  length, 

X  is  the  optical  wavelength,  and 

a  2  =  log  amplitude  variance  of  the 
signal. 

These  equations  may  be  related  to  the  log  intensity  variance 
by  noting  that 

V  *  4  V- 

We  then  have  the  following: 

Plane  Waves  a  2  =  1,24  K7/6  C  2  L1]/s 

z  n  . 

Spherical  Waves  a  2  =  0.52  K7/6  C  2  L11/6  . 
r  in 

The  restriction  that  ZQ  <<  (XL)  equates  to  a  minimum 

vortex  size  of  approximately  0.8  cm  for  X  =  .6328  micrometer 
(HeNe  laser  source)  with  a  path  length  of  100  meters.  More 
recent  observations  [Ref,  5]  report  general  agreement  with 
these  equations  for 

<?x 2  <  0.3 

and  a  2  <  1.2. 
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This  is  the  criteria  we  chose  to  use  in  defining  the  weak 
turbulence  regime  for  our  experiment. 

Direct  path  optical  scintillation  is  an  important 
consideration  in  many  practical  situations  of  military  interest 
Optical  laser  communication,  laser  spotting,  target  illumina¬ 
tion,  and  power  transmission  are  just  a  few  of  these. 

Rigorous  theoretical  analyses  of  direct  path  scintillation 
by  many  authors  over  the  last  few  decades  have  resulted  in 
a  fairly  high  degree  of  confidence  that  the  factors  involved 
in  direct  path  scintillation,  especially  turbulence,  are 
well  understood.  Theoretical  predictions  of  signal  log 
intensity  variance,  or  scintillation  strength,  agree  reasonably 
well  with  experimental  findings.  The  correlation  between 
direct  and  folded  path  scintillation  strength  is  less  well 
documented  or  understood.  In  his  doctoral  thesis  at  the 
Naval  Postgraduate  School,  Dr.  Avihu  Ze'evi  examined  the 
special  case  of  folded  path  scintillation  using  a  plane 
mirror  target  and  a  point  detector  in  the  vicinity  of  or 
at  the  same  location  as  the  optical  transmitter.  The 
correlation  of  the  scintillation  strength  for  a  folded  path 
with  the  scintillation  strength  for  a  direct  path  was 
formulated.  It  is  this  correlation  that  we  sought  to 
verify  in  our  experiment. 

B.  FOLDED  PATH  THEORY 

In  his  theoretical  treatment  of  the  folded  path  problem, 


1.  The  end  points  of  the  folded  path  are  close  in 
comparison  with  the  length  of  the  path. 
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2.  The  target  that  folds  the  path  is  a  specular  reflector 
such  that  there  is  no  discontinuity  of  the  light  rays. 

3.  An  infinite  plane  mirror  or  at  least  a  very  large  mirror 
is  used  as  tne  specular  reflector  so  that  all  of  the 
beam  is  reflected.  This  assumption  excludes  the 
effect  of  beam  wander  and  neglects  diffraction  effects. 

He  then  goes  on  to  state  that  the  basic  question  that  must 

be  addressed  in  determining  any  correlation  between  the 

direct  path  and  folded  path  situations  is  whether  or  not 

the  two  parts  of  the  folded  path  are  statistically  independent. 

If  one  can  assume  statistical  independence,  then  the  two 

parts  of  the  folded  path  can  be  treated  separately,  applying 

the  well  known  results  for  direct  path  propagation  to  each 

path.  Figure  1  illustrates  the  situation  where  the  light 

from  a  laser  source  is  partially  reflected  and  partially 

transmitted  at  a  specular  target.  <Sp2>  and  <S2L2:-  represent 

the  variance  of  the  signal  log  intensity  at  a  point  along 

the  folded  and  the  direct  paths,  respectively,  such  that 

the  distances  to  each  point  along  these  respective  paths 

are  equal.  Then  for  the  two  parts  of  the  folded  path  to 

be  statistically  independent, 

<sp  2>  -  <s2l2>. 

In  his  work,  Ze'evi  shows  that 

<Sp2  >  =  <S2l2>  Cl  +  f(yQ,  3)) 
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where  f(y0»  0)  is  the  path  dependency  function,  arrived  at 
by  numerical  integration.  The  path  dependency  function  is 
illustrated  in  Figure  2.  The  parameter,  y  ,  is  defined  to 
be 

,  L  Sin  <29) 

y° 

2  °/Lo» 

where  D  is  the  source  detector  displacement  and  0  is  the 
half  angle  between  the  two  parts  of  the  folded  path, 
illustrated  in  Figure  3.  For  small  angles  of  0,  the  f 
dependency  is  relatively  weak,  and  the  expression  for  <Sp2> 
then  becomes 

<Sp2>  =  <S2l2>  f1  + 

We  then  have  the  following: 


lim  f (y0,9)  *  1 


when  0  - -  0 ,  y - »  0 

J  o 


lim  f(yo,0)  »  0 
0  - *  IT 

T 


then 


<sf2(l,o)> 


2<s2l2>,  0=0 
<S2L2>*  0  =  I 


Thus,  it  is  clear  that  the  paths  are  statistically  independent 
only  for  2 0  =  however,  because  f(yQ,0)  decreases  rapidly 


Field  Source 


gure  3.  Diagram  of  the  System  Geometry  Showing  the 
Relationship  of  the  Quantities  D,  R,  and  29 
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for  small  values  of  yQ  and  approaches  slowly  to  zero,  Ze’evi 
defines  the  region  in  which  statistical  independence  exists 
as  that  where  f(yo,0)  =  e’1.  This  criteria  is  met  for  small 
9  where 

y  >4 
7  o  — 


and  statistical  independence  of  the  folded  path  now  depends 
solely  on  the  source  detector  displacement,  D. 

The  next  step  is  to  relate  the  displacement  to  the 
variances  of  the  folded  path  and  the  direct  path,  <Sp2>  and 
<Sl2>,  respectively.  We  are  not  interested  in  the  case  of 
<S2^2>  in  most  situations  of  practical  application.  It  is 
convenient  at  this  point  to  remember  that  the  variances  we 
are  working  with  are  generally  accepted  to  be  taken  from 
log  normal  distributions  in  the  weak  turbulence  regime,  and 
so  we  may  adopt  the  notation  a£2(L)  and  cj£2(F)  for  the 
direct  path  and  folded  path  variances,  respectively.  The 
subscript  l  indicates  a  log  normal  distribution. 

Using  the  results  for  the  log  intensity  variance  of 
spherical  waves  calculated  by  Tatarski  [Refs.  3,  4] 
as  follows: 


and 


where 


a£2  (L)  =  0.52  K  •?/s  Cn2  L1^ 


<S  I>  ' 


Ai^  r  (5/6)  L_ 

2  *3  L0j 


=  .  053  2  ^3  -r  3/2  K2  Cn2  L09/3 


r  (11/16) 
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Ze'evi  arrives  at  the  following  expression  for  a  path  of 
length  2L  in  the  weak  turbulence  regime: 

ot2(2L)  -  Aj  2s/6  Ls/6  <Szl2>. 

Relating  this  to  the  expression 

<SF2>  =  <S2L2>C1  +  f(yo» 

Ze’evi  shows  that 

o^OO  3  \  2^  L5/6  <Sl2>  Cl  +  f(yQ,  9)) 

and  finally, 

2 

R  =  ii  =  21*  (1  *  f(y  e). 

R  is  the  ratio  of  the  folded  path  scintillation  strength  to 
the  direct  path  scintillation  strength  in  the  weak  turbulence 
regime  (a^2  <  1.  2),  and  is  the  parameter  we  attempted  to 
verify  during  our  experiment.  The  range  of  R  is 

Z1^8  <  R  £  21*6 

or 

3.56  <_  R  <_  7 . 12  , 

By  restricting  ourselves  to  the  weak  turbulence  regime 

in  these  discussions,  we  have  ignored  such  effects  as 
saturation  of  the  scintillation  strength  and  the  transverse 


correlation  length  (the  smallest  lateral  distance  at  the 
receiver  that  would  allow  independent  intensity  fluctuations 
to  be  recorded)  of  the  signal.  In  addition,  because  we 
intended  to  measure  the  scintillation  strength  from  an 
(approximate)  point  detector  in  our  experiment,  we  have 
also  omitted  any  reference  to  aperture  averaging  of  the 
received  signal  intensity.  This  is  appropriate  in  keeping 
with  the  conditions  under  which  Ze'evi’s  folded  path  theory 
was  formulated.  With  these  caveats  stated,  we  may  proceed 
to  a  discussion  of  our  experimental  approach. 
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II.  EXPERIMENTAL  APPROACH 


A.  PROBLEM  DESCRIPTION 

An  initial  attempt  at  verifying  Ze'evi’s  theory  was 
attempted  by  Hodgini  [Ref.  6]  in  the  early  part  of  1982. 

This  was  of  benefit  to  us  in  our  work  because  much  of  the 
equipment  necessary  to  perform  the  experiment  was  available, 
and  because  several  factors  involved  in  the  experiment, 
e.g.,  the  stability  and  correlation  of  the  helium-neon  sources 
had  already  been  investigated. 

In  order  to  compare  the  correlation  between  the  direct 
path  and  folded  path  variance,  it  was  necessary  to  construct 
an  electro-optical  data  collection  system  that  could  measure 
the  desired  optical  parameters  (the  intensity  fluctuations 
from  the  two  different  path  geometries)  and  then  store  the 
data  in  a  histogram  type  of  arrangement  until  the  statistics 
of  the  intensity  distributions  could  be  calculated  and 
recorded.  The  primary  requirement  for  the  experiment  in 
order  that  the  direct  path  and  folded  path  scintillation 
strength  could  be  compared,  was  that  the  two  lasers  propagate 
over  the  same  uniform  surface  leading  to  the  detector.  For 
the  two  sources  to  be  subject  to  the  same  atmospheric 
conditions  over  the  path  length  meant  that  the  lasers  would 
have  to  be  alternately  pulsed  at  a  rate  high  enough  that 
the  turbulence  along  the  path  was  frozen  during  the  length 
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of  time  necessary  to  transmit  and  receive  both  pulses.  The 
atmosphere  can  be  considered  frozen  for  time  intervals  of 
up  to  1  millisecond  [Ref.  1].  A  basic  cycle  of  1000  Hertz 
was  chosen  giving  a  one  millisecond  period  wherein  first  the 
direct  then  the  folded  path  lasers  would  be  pulsed  and  the 
signals  sampled.  This  allowed  a  separation  of  approximately 
one  quarter  millisecond  between  each  successive  pulse. 

The  optical  setup  we  devised  is  illustrated  in  Figure  4, 
where  both  the  zero  displacement  case  and  an  arbitrary 
displacement  example  are  shown.  The  target  end  consists  of 
the  direct  path  laser,  whose  beam  is  modulated  by  a  mechanical 
chopper,  expanded  by  a  lens  system,  and  reflected  to  the 
detector  aperture  by  a  beam  splitter.  Behind  the  beam 
splitter  is  the  mirror  for  the  folded  path.  The  detector 
end  includes  both  the  folded  path  laser  with  an  optical 
modulation  system  (Pockels  cell  modulator  and  polarizer) , 
another  beam  splitter  for  the  case  where  the  folded  and 
direct  paths  are  directly  overlayed,  and  the  detector, 
dewar,  and  aperture.  The  signal  pulses  from  the  laser 
at  the  detector  end  then  propagate  down  range  to  the  target 
end,  pass  through  the  splitter,  reflect  off  of  the  mirror, 
and  are  returned  to  the  detector  via  the  same  path  as  the 
direct  path  laser. 

Our  two  greatest  problems  in  developing  the  experimental 
setup  were  timing  and  alignment.  Timing  was  of  the  utmost 
importance  due  to  the  rapidity  at  which  the  data  was  sampled 


23 


Direct  Path 


Figure  4.  Equipment  Diagram  for  D  =  Various  Displacements  (Top) 


and  because  of  the  need  to  store  the  sampled  data  alternately 
at  different  locations  corresponding  to  the  different 
intensity  distributions  from  the  two  paths.  Alignment  was 
likewise  critical  because  of  the  requirement  that  the 
optical  pulses  from  the  two  separate  sources  traverse  the 
same  physical  path  to  the  aperture  of  the  detector. 

In  consideration  of  these  two  difficult  requirements  and 
the  inherent  complexity  of  a  large  number  of  electronic  and 
optical  components,  the  system  was  first  assembled  and 
tested  in  the  basement  of  Spanagel  Hall.  A  laser  range 
running  along  the  ceiling  provided  a  145  meter  path  with 
eight  large  ventillator  ducts  delivering  heated  air  into 
the  optical  path  to  provide  the  necessary  turbulence. 

Once  system  development  was  complete,  a  variety  of 
interior  data  runs  were  conducted  for  final  system  checkout. 
The  entire  system,  exactly  as  constructed  and  tested  inside, 
was  then  relocated  to  an  outside  site.  This  site  was  a  100 
meter  long,  flat  grassy  area  alongside  Spanagel  Hall  and 
approximately  30  meters  away.  All  optical  components  were 
situated  on  concrete  tables  and  anchored  with  lead  bricks 
for  stability,  while  the  signal  processing  equipment  and 
computers  were  placed  on  movable  tables.  Later,  another 
outdoor  site  atop  Spanagel  Hall  would  be  tried  that  provided 
a  shorter  path  but  more  scintillation. 


B.  OPTICAL  EQUIPMENT 


The  final  system  consisted  of  one  optical  configuration 
for  all  data  runs  except  the  zero  displacement  run,  which 
incorporated  an  additional  beam  splitter  at  the  detector 
end.  Both  of  these  configurations  are  depicted  in  Figure  4. 

1 .  Laser  Sources 

The  two  lasers  utilized  throughout  the  experiment 
were  Spectra  Physics  HeNe  lasers  with  a  wavelength  of  .6328 
micrometer  and  an  output  of  .95  milliwatt.  These  lasers 
were  chosen  as  sources  for  the  following  reasons:  One,  they 
produce  light  in  the  visible  range,  making  alignment  easier; 
two,  the  output  power  and  therefore  the  hazard  level  was 
relatively  low;  and  three,  the  correlation  between  these  two 
lasers  when  pulsed  alternately  with  a  pulse  separation  of 
0.3  millisecond  over  the  same  optical  path,  through  the  same 
turbulence  was  excellent.  Hodgini  reported  source  variances 
for  these  lasers  as  less  than  3xl0-4,  and  the  correlation 
coefficient  (ratio  of  source  variances)  for  the  two  HeNe 
lasers  as  .996  [Ref.  6],  Although  the  lasers  had  a  tendency 
to  become  unstable  with  time,  the  correlation  between  the 
two  sources  did  not,  and  therefore  they  could  be  used 
together  in  measuring  scintillation  with  some  confidence 
that  their  internal  stability  (source  variance)  would  not 
be  a  contributor  to  the  scintillation  we  were  measuring. 


2.  Target  End  Optics 

The  target  end  optics  consisted  of  a  beam  expanding 
lens  system,  a  set  of  adjustable  polarizers  for  fine  intensity 
adjustment  of  the  direct  path  laser  signal  intensity,  a  four 
inch  optical  beam  splitter,  and  mirror.  The  lens  system  was 
required  in  order  to  expand  the  direct  path  laser  spot  size 
to  the  approximate  size  of  the  folded  path  laser  spot  at 
the  target  end  beam  splitter.  This  also  helped  match  folded 
and  direct  path  signal  intensities  at  the  detector  aperture, 
and  made  system  alignment  easier.  Both  the  direct  path 
signal  and  the  folded  path  signal  would  then  transit  through 
the  same  volume  of  atmosphere  in  reaching  the  detector 
aperture.  The  first  lens,  actually  a  microscope  objective, 
received  the  laser  beam  and  focused  it  to  a  point  near  the 
focal  point  of  the  second  lens,  which  then  collimated  the 
beam  entering  the  beam  splitter.  Two  polarizers,  placed 
ahead  of  the  microscope  objective,  were  used  to  make  final 
reductions  in  the  direct  path  signal  intensity  to  closely 
match  the  intensity  of  the  folded  path  signal. 

The  four  inch  beam  splitter  was  a  critical  component 
because  the  folded  path  signal  had  to  pass  through  it  twice 
and  the  direct  path  signal  had  to  be  reflected  by  it  enroute 
to  the  detector.  It  was  therefore  of  extreme  importance 
that  no  structure  be  imparted  to  the  transiting  pulses  and 
also  that  reflection  losses  be  minimal  so  that  the  signal 
intensities  reaching  the  detector  be  as  high  as  possible. 


To  satisfy  these  conditions,  a  four  inch  diameter,  half  inch 
thick  optical  window  was  acquired  with  an  advertised  flat¬ 
ness  correct  to  1/10  wavelength  and  parallelism  within  two 
seconds  of  arc.  An  antiref lective  coating  was  applied  and 
one  side  of  the  splitter  silvered  to  provide  50  percent 
transmission,  the  latter  meaning  50  percent  of  the  signal 
intensity  from  either  source  was  lost  at  each  pass  through 
the  beam  splitter.  Coating  and  silvering  the  splitter  also 
produced  the  beneficial  result  that  the  maximum  signal 
intensity  from  the  direct  path  laser  reflected  by  the  back 
side  of  the  beam  splitter  was  less  than  0.4  percent  of  the 
signal  intensity  reflected  from  the  front  of  the  splitter, 
and  therefore  could  be  neglected.  The  target  mirror  utilized 
was  a  five  inch  optical  flat.  The  important  qualities  required 
were  good  optical  flatness  and  large  size.  The  flatness 
had  been  previously  verified  using  an  interferometer  and 
the  five  inch  diameter  proved  large  enough  to  accomodate 
the  folded  path  spot  size  after  its  initial  100  meter  transit. 

C.  PULSE  FORMING  AND  DETECTION  EQUIPMENT 

The  pulse  forming  equipment  in  our  experiment  included  a 
mechanical  chopper  to  modulate  the  direct  path  laser  and 
and  electro -optical  modulation  system  to  modulate  the  folded 
path  laser.  After  the  signals  arrived  at  the  aperture,  they 
were  converted  to  an  analog  voltage  signal  by  the  detector 
and  amplified  before  being  log  converted  and  stored. 
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1.  Mechanical  Chopper 

The  mechanical  chopper  was  positioned  between  the 
second  polarizer  and  the  microscope  objective  at  the  target 
end.  A  circular  wheel  with  an  open  to  closed  ratio  of  1:4 
was  rotated  using  a  3000  rpm  electrical  motor,  resulting  in 
a  chopping  frequency  of  1000  Hertz.  A  light  emitting  diode 
mounted  opposite  a  detector  was  also  placed  on  the  chopper 
and  provided  the  reference  trigger  signal  for  the  entire 
experimental  system.  This  trigger  was  led  to  an  Interstate 
P12  pulse  generator  which  delayed  the  signal  to  provide  the 
trigger  for  the  modulation  system  on  the  folded  path  laser. 

2 .  Modulation  System 

A  Coherent  Associates  Model*  3003  modulation  system 
(Pockels  cell)  was  used  for  pulse  forming  of  the  folded 
path  laser.  After  the  pulse  of  the  direct  path  laser  was 
initiated  by  the  mechanical  chopper,  a  delay  of  approximately 
.4  millisecond  was  provided  by  the  P12  described  above, 
and  this  delayed  trigger  was  led  to  the  Pockels  cell  modula¬ 
tion  system  to  initiate  the  folded  path  pulse.  Folded  path 
pulse  width  could  be  matched  to  that  provided  by  the  mechanical 
chopper  by  varying  the  pulse  width  at  the  P12. 

3 .  Detector 

The  detector  utilized  for  our  intensity  measurements 
was  a  Model  50  EHS  Silicon  Avalanche  Detector  manufactured 


by  the  General  Electric  Company.  When  operated  at  voltages 
of  between  1200  and  1500  volts  and  cooled  in  a  dewar  of 


liquid  nitrogen,  the  detector  provided  an  extremely  high 
signal  to  noise  ratio.  The  high  voltage  caused  high 
amplification  of  weak  signals  by  placing  the  detector  into 
its  avalanche  region,  and  the  cooling  suppressed  noise.  A 
high  signal  to  noise  ratio  was  needed  because  of  the  extremely 
low  signal  power  available  at  the  detector  aperture  from  the 
HeNe  sources.  Although  we  did  not  verify  the  linearity  of 
this  detector  over  the  dynamic  range  of  incident  signal 
intensity  available  to  us,  the  detector's  linear  response  in 
the  weak  signal  regime  had  been  previously  established  by 
the  atmospheric  physics  group  at  NPS  and  we  neglected  any 
nonlinearity  at  this  point. 

In  order  to  make  our  detector  system  a  point  detector, 
the  aperture  needed  to  be  as  small  as  possible  and  still 
pass  enough  light  to  enable  the  detector  to  register  a 
signal.  Hodgini  [Ref.  6]  shows  that  the  Fresnel  zone  site 
at  the  aperture  is  the  most  important  criteria.  Zone  size 
can  be  determined  as  follows: 

=  (XL)1/2 

where  is  the  Fresnel  zone  radius, 

r 

X  =  .6328  micrometer,  and 


L  is  the  range. 
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The  Fresnel  zone  size  at  145  meters  is  9.58  mm,  at  100 
meters  is  7.95  mm,  and  at  40  meters  is  5.03  ram.  This  is  a 
smaller  length  than  the  lateral  coherence  length  of  the 
light,  either  outside  or  in  the  lab,  and  therefore  an 
aperture  size  smaller  than  the  Fresnel  zone  size  for  a 
particular  range  would  approximate  a  point  detector.  Using 
this  constraint,  Hodgini  had  settled  on  a  3  mm  diameter 
aperture  and  we  adopted  this  size  as  well.  Light  passing 
through  this  opening  was  then  focused  on  the  detector  by  a 
single  convergent  lens  placed  a  focal  length  from  the 
detector  surface.  We  did  not  use  any  spectral  filters 
between  the  aperture  and  the  lens  as  all  of  our  experimental 
work  was  done  under  conditions  of  very  low  ambient  light. 

4.  Signal  Amplifier 

The  output  from  the  detector  was  led  directly  to  a 
Model  113  Princeton  Applied  Research  (PAR)  Low  Noise  Amplifier 
This  device  provided  high  signal  gain  and  an  adjustable 
bandwidth  with  high  and  low  frequency  roll-offs.  These 
features  were  necessary  to  scale  the  amplitude  of  the  signal 
reaching  the  analog  to  digital  converter  on  the  storage 
computer,  and  also  to  obtain  the  optimum  signal  shape.  One 
point  in  particular  should  be  made  regarding  the  selection 
of  signal  bandwidth,  and  that  is  this:  we  had  to  be  very 
careful  that  the  low  frequency  roll-off  was  not  set  too 
high,  resulting  in  a  dip  in  the  background  level  following 


the  signal  (removal  of  too  many  Fourier  components  from  the 
wave  shape)  and  leading  to  an  erroneous  background  sample  in 
the  signal  processing  performed  by  the  demodulators.  In 
other  words,  we  wanted  to  eliminate  as  much  noise  as  possible 
but  not  damage  the  optimum  signal  shape  in  doing  so. 

D.  SIGNAL  PROCESSING  EQUIPMENT 

Scintillation  effects  from  the  atmosphere  along  the  path 
traveled  by  the  light  pulses  caused  the  intensity  at  the 
aperture  and  therefore  at  the  detector  to  be  a  function  of 
time.  This  same  effect  was  also  present  in  the  background 
intensity  just  before  and  after  the  signal.  We  wanted  to 
be  able  to  sample  the  largest  instantaneous  value  of  the 
intensity  during  the  pulse,  and  the  value  of  the  background 
intensity  following  the  pulse.  The  difference  between 
these  two  values  would  be  the  true  signal  intensity.  This 
process  of  sampling  and  returning  the  absolute  intensity  is 
termed  "demodulation",  and  the  equipment  we  needed  to  perform 
this  task  for  both  laser  signals  was  available  in  the  form 
of  two  demodulators  used  by  the  atmospheric  physics  group 
at  NPS.  Additionally,  since  the  actual  measure  of  the 
scintillation  strength  is  the  log  intensity  variance,  the 
absolute  intensity  magnitudes  had  to  be  converted  to  logarithm 
values  prior  to  being  digitized  and  stored.  Our  signal 
processing  equipment  then  consisted  of  two  demodulators, 
two  log  converters,  and  the  electronic  switch. 


1.  Demodulators 


Two  demodulators  were  utilized  in  our  experiment, 
one  for  both  the  direct  path  and  the  folded  path  signal. 

The  function  of  each  was  identical.  A  demodulator  received 
the  signal  directly  from  the  PAR  low  noise  amplifier  described 
previously.  A  trigger  pulse  activated  the  demodulator  to 
receive  only  its  respective  input  signal,  after  which  the 
demodulator  had  to  accomplish  four  actions  to  output  a 
differential  signal.  Figure  5  illustrates  these  actions 
for  the  direct  path  signal;  the  folded  path  signal  was 
exactly  the  same,  simply  occurring  about  0.3  millisecond 
later.  The  stretch  cutoff  function  took  the  instantaneous 
maximum  of  an  incoming  signal  and  held  it  for  an  adjustable 
period  of  time,  allowing  the  signal  sample  to  be  set  within 
this  time  frame.  After  the  stretch  cutoff  was  released, 
the  demodulator  saw  only  the  background  signal,  and  a 
background  sample  was  taken.  The  demodulator  electronically 
took  the  difference  between  the  two  values  it  had  sampled 
and  output  the  difference.  This  differential  signal  was 
held  until  the  next  input  and  input  trigger  arrived. 

2 .  Log  Converters 

The  differential  output  of  each  demodulator  was 
sent  to  its  own  Hewlett-Packard  7562A  Log  Converter.  These 
devices  produced  a  dc  output  voltage  that  was  proportional 
to  the  log  of  the  dc  input  voltage  in  the  following  manner; 
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where  Vpis  the  demodulator  differential  output,  VL  is  the 
log  output,  and  a  and  b  were  constants  determined  by  the 
design  of  the  equipment.  The  constant,  a,  simply  resulted 
in  a  shift  of  the  mean  of  the  distribution,  and  could  be 
different  for  each  log  converter.  This  was  identical  to 
a  uniform  increase  in  incident  signal  intensity  and  did  not 
affect  the  log  variance,  as  shown  in  Chapter  I.  The 
constant,  b,  had  to  be  the  same  for  each  log  converter, 
or  the  variance  ratio 

•('W 

would  be  multiplied  by  a  gain  factor.  We  measured  b  for 
each  log  converter  and  found  the  values  to  be  within  .009 
of  each  other.  This  confirmed  the  accuracy  of  the  two 
parallel  data  lines  reaching  the  electronic  switch.  The 
two  log  converted  differential  signals  corresponding  to 
the  direct  path  and  folded  path  signal  intensities  were 
now  ready  for  digitizing  and  storage  in  the  data  storage 
computer . 

3 .  Electronic  Signal  Switching 

The  Nicolet  1080  Stored  Program  Computer  was  used 


for  data  storage  during  an  experimental  run.  Each  experimental 
run  consisted  of  collecting  sixteen  pairs  of  variance  values 
at  a  specific  value  of  the  displacement,  D,  and  between  10 


to  13  displacements  comprised  an  experimental  try  at 
validating  the  folded  path  theory.  Because  the  NIC- 80  had 
only  one  analog  to  digital  converter  (A DC)  port,  we  time- 
shared  the  port  between  the  two  input  signal  trains.  A 
new  signal  arrived  at  the  ADC  every  0.3  millisecond.  A 
switch  was  required  that  would  receive  the  output  signals 
from  the  log  converters  and  sequence  these  signals  to  the 
ADC.  Our  switch  was  built  using  the  AD7510D1  analog  chip. 
The  function  of  the  electronic  switch  is  illustrated 
graphically  in  Figure  6,  where  the  two  independent  data 
signals,  represented  by  dotted  lines,  are  shown  being 
sequenced  through  the  electronic  switch  to  the  NIC-80.  The 
solid  lines  represent  system  timing  signals.  An  Interstate 
Electronics  Corporation  P25  Pulse  Generator  was  synchronized 
to  the  delayed  trigger  of  the  folded  path  optical  modulation 
system  and  output  two  continuous  pulse  trains,  one  the 
complement  of  the  other.  Both  the  positive  and  negative 
pulse  trains  went  to  the  electronic  switch  and  controlled 
the  transmission  of  the  alternating  signals  to  the  ADC.  In 
addition,  the  positive  pulse  train  went  to  the  port  on  the 
NIC-80  labeled  "Sense  2"  to  tell  the  computer  in  which  data 
storage  area  to  assign  the  incoming  data  point.  When  a 
positive  value  was  present  at  the  switch  and  Sense  2,  the 
folded  path  log  intensity  signal  was  allowed  to  transit 
the  switch,  be  digitized  in  the  ADC,  and  enter  the  folded 


path  storage  area  in  the  computer.  The  relationship  of 
these  signals  in  time  is  shown  in  Figure  7. 

E.  COMPUTERS 

The  computer  system  which  we  needed  for  our  experiment 
was  available  and  already  adapted  by  the  atmospheric  physics 
group  at  NPS  for  the  express  purpose  of  storing  and  processing 
the  high  speed  analog  signals  from  direct  path  scintillation 
measurements.  This  system  consisted  of  three  components, 
the  Nicolet  Instrument  Corporation  Model  NIC-80,  the  Hewlett- 
Packard  (HP)  9825A  Calculator,  and  the  HP  9871A  Printer. 

To  these  three  we  added  a  fourth  component,  the  HP  9862A 
Calculator  Plotter,  to  obtain  a  graphical  output  of  the 
data  as  an  aid  in  determining  its  validity.  Each  component 
of  the  system  is  discussed  in  detail  in  the  paragraphs  that 
follow. 

1.  NIC-80 

The  NIC-80  is  a  general  purpose  computer,  highly 
suited  to  the  role  of  high  speed  data  storage  and  manipulation. 
It  is  configured  for  use  with  a  CRT  and  has  a  hardwired 
analog  to  digital  converter  (ADC)  as  well  as  several  ports 
for  the  purpose  of  timing  and  control.  Signals  from  both 
laser  sources  were  sent  to  the  NIC-80  via  the  electronic 
switch,  arriving  at  the  ADC  where  they  were  converted  to  a 
digital  value.  This  value  was  scaled  by  the  stored  program 
to  correspond  to  a  range  of  memory  addresses  in  the  machine, 
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and  each  time  the  value  occurred,  a  count  was  added  at  the 
corresponding  memory  address.  Only  positive  values  from  the 
ADC  were  counted  by  the  stored  program,  which  helped  smooth 
out  the  tail  regions  of  each  distribution.  During  a  data 
collection  period,  each  source  had  its  own  block  of  memory 
addresses  and  thus  a  histogram  of  counts  versus  point  signal 
intensity  at  the  detector  was  developed  for  each  laser  source. 
The  range  of  each  histogram  (intensity  distribution)  in  the 
NIC-80  was  256  memory  spaces  or  "bins",  corresponding  to  a 
range  of  0  to  256  from  the  ADC.  Through  tests  run  on  artificial 
data,  it  was  determined  that  the  bin  widths  at  the  ADC,  that  is, 
the  analog  signal  range  per  value  (or  bin)  at  the  ADC  was 
so  close  to  a  constant  value,  that  for  the  purposes  of  our 
statistical  analysis  of  each  source’s  histogram,  we  could 
assume  equal  bin  widths  of  signal  intensity.  This  proved 
to  be  an  excellent  assumption,  and  enormously  simplified  the 
statistical  calculations  that  followed.  The  smoothness  of 
the  intensity  distribution  that  we  got  from  this  process 
of  digitizing  and  storage  was  almost  entirely  a  function 
of  the  average  bin  width  in  the  ADC.  Our  average  bin  width 
was  33.169  millivolts,  and  this  proved  adequate  for  the  range 
of  signal  intensities  that  we  used.  The  value  of  the 
original  signal  intensity  at  the  detector  was,  of  course, 
not  important,  simply  the  variance  of  mean  signal  intensity. 

Each  run  consisted  of  16,384  samples  per  laser  source  or 
32,768  samples  total,  and  took  approximately  15  seconds  to 


complete.  After  completing  a  data  run,  the  NIC-80  program 
would  reassign  a  new  block  of  memory  for  each  laser  source 
and  the  process  would  be  repeated  until  16  runs  or  32  data 
pairs  had  been  recorded.  The  entire  16  run  period  took 
approximately  four  minutes,  during  which  time  it  was  possible 
to  monitor  visually  on  the  CRT  the  building  up  of  the 
histogram  for  each  source.  In  other  words,  we  had  a  real 
time  picture  of  the  scintillation  data  during  each  run. 

This  was  very  useful  for  monitoring  the  relative  variances 
of  each  source,  as  well  as  their  relative  mean  intensities, 
with  the  result  that  the  quality  of  the  data  was  immediately 
apparent  and  could  be  assessed  before  devoting  too  much  time 
to  a  run.  Such  factors  as  poor  alignment  and  faulty  system 
timing  could  be  spotted  immediately  from  the  nature  of  the 
intensity  distributions  displayed  on  the  CRT. 

The  NIC-80  was  programmed  entirely  in  assembly 
language,  compiled  and  loaded  from  the  HP-9825A.  The 
assembly  language  program  for  the  NIC-80  is  included  in 
Appendix  A.  Timing  pulses  arriving  at  the  ports  labeled 
Sense  1  and  Sense  2  were  interpreted  by  the  computer  program 
to  tell  the  machine  which  source's  signal  was  currently 
present  at  the  ADC  and  thus  the  program  directed  the  storage 
of  the  data  count  at  the  proper  memory  address.  For  example, 
a  high  voltage  present  at  Sense  2  meant  that  the  signal  from 
the  direct  path  source  was  ready  for  digitizing,  and  a  low 
voltage  meant  that  the  signal  was  from  the  folded  path  source. 
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A  concurrent  trigger  pulse  at  Sense  1  then  allowed  the  signal 
at  the  ADC  to  be  digitized.  Sense  1  therefore  had  to  be 
triggered  at  twice  the  frequency  of  the  data  rate  from  either 
source.  At  first  glance  this  seems  redundant  with  the 
function  of  Sense  2,  but  it  really  is  not;  Sense  1  paced  the 
digitizing  and  Sense  2  determined  which  laser  source  was 
being  digitized. 

2 .  HP-9825A  Calculator 

The  HP-9825A  Calculator  was  the  operator-machine 
interface  for  our  experiment,  and  the  program  for  the  9825A 
controlled  the  function  of  the  NIC-80,  the  printer,  and  the 
plotter.  Control  of  the  NIC-80  was  accomplished  through 
a  computer  interface  built  by  Processor  Edmund  Milne  at  NPS , 
and  the  control  signals  consisted  of  pairs  of  binary  values 
exchanged  as  initiation  and  answer  "handshakes"  at  critical 
points  to  control  the  execution  of  the  NIC-80  program.  The 
9825A  program  was  written  in  HPL  and  was  initially  structured 
to  collect  all  the  data  for  a  given  value  of  D  in  the  NIC-80 
then  transfer  the  data  to  the  9825A  for  processing  and 
graphing  immediately  after  the  run.  Approximately  30  minutes 
was  then  required  to  print  and  graph  the  data  for  all  16  pairs 
this  was  the  method  used  to  collect  the  data  for  -  .perimental 
results  one  and  two,  discussed  in  Chapter  III.  A  copy  of 
this  program  for  the  9825A  appears  in  Appendix  B.  For  reasons 
to  be  discussed  later,  we  decided  late  in  the  project  to 
try  a  different  method  of  data  collection  and  processing 


in  which  it  was  desired  to  shorten  the  elapsed  time  between 
runs,  time  which  previously  had  been  devoted  to  printing  and 
graphing.  Our  goal  was  to  try  and  collect  fewer  numbers  of 
data  pairs  at  each  value  of  D,  and  to  try  and  collect  all 
the  data  for  all  values  of  D  that  we  intended  to  use  --12 
or  13  points  usually--in  30  minutes  or  less  by  recording  the 
data  from  the  NIC-80  on  magnetic  tape  and  moving  to  the  next 
displacement  without  processing  the  data.  This  was  the 
method  used  to  record  the  data  from  our  third  experiment, 
and  the  program  to  do  this  is  also  included  in  Appendix  B. 

One  final  note  needs  to  be  made  regarding  the  9825A  interface 
with  the  NIC-80.  The  NIC-80  data  word  is  20  bits  long,  but 
the  parallel  input/output  bus  is  16  bits  long.  This  necessitated 

the  transfer  of  data  from  the  NIC-80  to  the  9825A  in  two  10 

bit  segments,  which  is  two  transfers  per  word  of  data.  The 
first  transfer  included  the  10  highest  bits  from  the  word 
and  the  second  transfer  included  the  10  lowest  bits  from 
the  word;  the  first  part  was  multiplied  by  1024  in  the  9825A 
and  added  to  the  second  part  to  form  the  complete  data  value. 

This  process  was  repeated  for  all  the  values  in  each  intensity 
distribution  from  both  laser  sources  in  a  given  run. 

After  loading  all  the  values  in  a  pair  of  intensity 

distributions  from  a  run,  the  9825A  was  now  ready  to 

calculate  the  statistical  parameters  that  we  desired.  This 
was  accomplished  in  a  straightforward  manner,  using  the 


accepted  formulas  for  mean,  variance,  standard  deviation,  and 
chi-square  shown  in  Appendix  C.  .  No  curve  fits  were  used 
and  no  data  points  were  excluded  to  reduce  the  effect  of 
noise  in  the  tails  of  the  log  normal  distributions.  Even 
so,  the  values  of  the  variance  calculated  by  our  program 
and  another  program  from  Professor  Edmund  Milne  that  used 
a  curve  fit  to  determine  the  variance  were  within  one  percent 
of  each  other.  This  strengthened  our  confidence  in  our 
approach.  In  order  to  determine  how  closely  the  log  normal 
curves  corresponded  to  a  true  Gaussian  curve,  we  chose  the 
chi-square  calculation  with  a  goodness  of  fit  parameter  of 
288.87  for  a  95%  confidence  interval,  derived  from  the  number 
of  degrees  of  freedom  (253)  in  our  data  histograms.  This 
portion  of  the  calculation  required  that  we  use  values  from 
the  Standard  Normal  Table,  and  we  interpolated  linearly 
between  values  in  the  table  out  to  the  sixth  decimal  place 
to  obtain  a  smooth,  continuous  true  normal  for  each  set  of 
data.  Test  runs  conducted  with  data  points  derived  from  the 
equation  of  a  Gaussian  curve  were  used  to  verify  the  statistical 
values  arrived  at  by  our  program,  including  the  chi-square 
value,  which  in  theory  would  be  zero  for  a  true  Gaussian 
distribution.  Different  runs  utilized  different  values  for 
the  mean,  variance,  and  maximum  ordinate;  in  all  cases  the 
program  determined  the  correct  values  and  returned  a  chi- 
square  value  less  than  0.3.  A  sample  test  run  is  shown  in 
Figure  8.  After  performing  this  software  validation,  we  had 

complete  confidence  in  the  programs. 
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Figure  8.  Chi-Square  Validation  Run 


The  HP - 9 8 7 1A  Printer 


The  HP-9871A  Printer  is  an  impact  printer  controlled 
by  the  9825A.  We  programmed  the  982 SA  to  print  automatically 
the  standard  deviation  and  chi-square  value  for  each  intensity 
distribution  in  a  run.  These  printouts  appear  in  Appendices 
D  through  F,  along  with  representative  plots  from  the 
runs.  A  header  on  each  page  lists  the  date  the  data  was 
recorded,  whether  or  not  the  detector  was  cooled  with  the 
aperture  in  place,  what  the  gain  was,  the  avalanche  voltage, 
bandpass,  and  the  value  of  the  displacement  in  cm  for  that 
particular  run.  In  the  columns  below  this  are  listed  the 
number  of  data  points  that  actually  made  it  into  the  histo¬ 
gram,  the  variance,  and  the  chi-square  value.  In  each  case, 
the  direct  path  parameters  are  displayed  first,  then  the 
folded  path  parameters. 

4.  HP-9862A  Calculator  Plotter 


The  last  component  in  our  computer  system  was  the 
HP-9862A  calculator  plotter.  We  used  the  plotter  to 
display  graphically  the  data  points  from  each  histogram, 
as  well  as  the  true  normal  curve  calculated  by  the  chi- 
square  approach.  The  smoothness  of  the  latter  curve  was 
a  result  of  the  linear  interpolation  of  the  normal  table 
referred  to  earlier.  Plotting  the  data  allowed  us  to  judge 
easily  the  quality  of  the  log  normal  assumption  for  the 
scintillation,  and  also  discover  where  the  data  differed 


from  the  true  normal  curve  if  the  chi-square  value  was  in 
excess  of  our  goodness  of  fit  parameter.  Plotting  the  data 
was  not  automatic  for  each  run  as  was  the  printout;  instead, 
we  took  one  representative  set  of  curves  for  each  value  of 
the  displacement  that  we  measured.  These  plots  are  arranged 
with  their  respective  printouts  in  Appendices  D  through  F. 
The  printout  for  a  run  overlays  the  two  plots  that  go  with 
it.  Note  that  we  did  not  always  take  the  best  set  of  data 
from  a  run  (the  lowest  chi-square  values),  but  instead  we 
chose  to  keep  plots  that  were  representative  of  the  entire 
run  as  a  whole.  Again,  each  plot  includes  a  header  coded 
as  in  the  printout. 


EXPERIMENTAL  WORK 


There  were  three  main  goals  of  our  experiment.  First, 
we  wanted  to  design  a  system  that  would  allow  us  to 
simultaneously  record  the  scintillation  over  both  a  direct 
and  a  folded  path.  Second,  we  wanted  to  calculate  the 
scintillation  strength  for  each  path,  which  meant  that  we 
needed  to  calculate  the  variance  of  the  received  intensity. 

In  addition,  we  wanted  to  know  after  we  had  the  variance  how 
accurate  had  been  our  assumption  of  a  log  normal  distribution 
Our  third  goal  was  to  compare  our  experimental  results  of 
the  relationship  between  the  direct  path  and  folded  path 
variances  to  the  theoretical  prediction  by  Ze'evi. 


B.  PREPARATION 

In  order  to  perfect  alignment  procedures  and  timing 
sequences,  a  large  number  of  runs  were  conducted  on  our 
inside  range  in  the  basement  of  Spanagel  Hall.  Due  to  the 
unknown  turbulence  conditions  in  the  basement  and  the 


differences  between  operating  in  the  open  atmosphere  and 
in  the  building,  these  runs  were  only  useful  for  checking 
out  the  optics,  timing,  and  data  reduction  techniques.  No 
data  from  these  runs  is  included  in  this  report.  We  should 
mention  at  this  point  that  we  did  not  undertake  to  measure 


the  overall  system’s  dynamic  range  or  linearity  before 
concluding  our  test  phase.  The  reason  for  this  is  that  all 
of  the  component  equipment  comprising  our  system,  including 
the  detector,  was  operated  well  within  the  linear  regime  for 
each.  We  thus  felt  that  further  investigation  of  these 
factors  was  neither  profitable  nor  necessary. 

C.  FIRST  OUTSIDE  EXPERIMENT 

We  transferred  our  entire  system,  including  the  target 
end  equipment,  detector  end  equipment,  signal  processing 
equipment,  and  the  computers,  outside  to  a  100  meter  range 
on  the  South  side  of  Spanagel  Hall  approximately  30  meters 
from  the  building.  We  arbitrarily  selected  100  meters  as 
the  range  because  we  felt  that  we  would  need  a  long  distance 
in  the  relatively  sheltered  area  south  of  the  building  in 
order  to  get  enough  scintillation  to  have  variances  above 
the  source  variance  regime  of  10"4.  We  went  out  30  meters 
from  the  base  of  the  building  to  try  to  get  away  from  any 
nonuniform  turbulence  effects  near  the  Hall.  It  quickly 
became  evident  that  daylight  runs  utilizing  the  HeNe  sources 
would  not  be  possible.  The  power  available  from  each  source 
was  not  enough  to  be  detectable  by  the  eye,  making  system 
alignment  almost  impossible.  Hodgini  notes  in  his  thesis 
that  optical  alignment  proved  to  be  of  extreme  importance 
in  obtaining  accurate  results,  and  for  this  reasons  we  wanted 
to  be  absolutely  certain  that  the  alignment  was  correct  and 


that  both  spots  from  the  laser  sources  were  centered  on 
the  detector  aperture  while  running  the  experiment.  Because 
of  the  small  size  of  the  aperture  (3mm)  and  the  presence  of 
the  optical  splitter  that  further  reduced  laser  signal 
intensity,  our  ability  to  align  the  system  visually  was 
critical;  for  this  reason  we  decided  to  do  the  experiment 
at  night. 

The  direct  path  laser  source  was  collimated  to  produce 
a  constant  spot  size.  This  was  done  in  order  to  match 
approximately  the  direct  and  folded  path  spot  sizes  at  the 
target  end  splitter  so  that  pulses  from  the  direct  path 
source  would  be  traveling  through  the  same  volume  of  atmosphe 
and  hence  subject  to  the  same  turbulence  conditions  as  the 
folded  path  pulses  on  their  return  to  the  detector.  The 
variance  in  the  log  normal  curve  resulting  from  the  direct 
path  source  would  then  be  larger  than  if  the  light  had  been 
allowed  to  diverge  normally  because  the  pulse  wave  fronts 
would  be  plane  instead  of  spherical.  We  corrected  for  this 
effect  when  analyzing  the  data. 

The  first  time  that  we  ran  the  experiment  outside,  11 
different  displacements  were  sampled  from  D  =  0  to  D  =  100  cm 
in  intervals  of  10  cm.  The  data  was  recorded  over  a  period 
of  eight  hours  the  night  of  3  November  between  8:30  PM  and 
4:30  AM.  Pictures  of  the  equipment  and  the  range  alongside 
Spanagel  Hall  appear  in  Figures  9  through  15.  The  weather 


Figure  9. 


An  Overview  of  the  Data  Reduction,  Signal  Processing, 
Folded  Path  Source  Equipment,  and  Detector  Assembly. 

The  leftmost  bay  contains  (from  top  to  bottom)  the 
HP-9825A  Calculator,  HP-9871A  Printer,  and  HP-9862A 
Calculator  Plotter.  To  the  right,  this  table  contains 
the  signal  processing  equipment  which  is  further 
identified  in  Figure  10.  The  optical  equipment  and 
detector  are  also  further  detailed  in  Figure  11.  The 
centermost  table  supports  the  PAR  amplifier  (top,  left) 
the  high  voltage  power  supply  for  the  detector  (top, 
right),  and  the  modulation  system  power  supply  (bottom) 


Figure  10. 

Left  bay,  top  to  bottom:  Direct  path  demodulator,  direct 
path  log  converter,  and  signal  generator. 

Middle  bay,  top  to  bottom:  Electronic  switch,  folded  path 
demodulator,  folded  path  log  converter,  two  P-12's  (leftmost 
sent  to  Sense  1,  rightmost  is  delayed  trigger  for  Pockels  cell, 
oscilloscope,  and  P-2S. 

Right  bay,  top  to  bottom:  HP-9825  to  Nicolet  1080  Interface, 
oscilloscope  for  Nicolet  1080  display,  and  Nicolet  1080. 
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Figure  11. 

Left  to  right  on  optical  bench:  detector  and  dewar  assembly 
with  3  mra  aperture  installed;  folded  path  Spectra  Physics 
HeNe  laser  with  pulse  modulation  system  (Pockels  Cell)  in 
front 
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Figure  12 


From  right  to  left:  direct  path  Spectra  Physics  laser, 
polarizer  pair,  mechanical  chopper,  microscope  objective, 
collimating  lens,  and  optical  splitter  with  target  mirror 
behind  it. 
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Figure  13.  Another  View  of  Direct  Path  End  Highlighting 
the  Optical  Splitter  and  Target  Mirror 
Arrangement 
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was  clear  and  cool  with  a  very  slight  wind,  and  no  problems 
were  encountered  in  system  alignment  or  operation.  The  log 
normal  distributions  for  each  run  appeared  very  reasonable, 
in  that  they  were  all  approximately  Gaussian,  with  a  chi- 
square  value  close  to  and  in  some  cases  less  than  the  required 
goodness  of  fit  value  of  288.87.  Note  that  though  the 
standard  deviations  of  the  direct  path  sources  were  usually 
greater  than  those  of  the  folded  path  source,  it  must  be 
remembered  that  the  direct  path  data  was  collected  from 
plane  wave  pulses,  and  the  folded  path  data  was  collected 
from  spherical  wave  pulses.  The  relationship  of  plane  wave 
scintillation  variance  to  spherical  wave  scintillation 
variance  can  be  deduced  from  the  ratio  of  the  equations  by 
Tatarski  in  Chapter  I,  and  is  0.52/1.24  or  .4194.  In  order 
to  take  the  ratio  of  the  folded  path  variance  to  the  direct 
path  variance,  we  first  converted  the  direct  path  variances 
by  multiplying  each  value  by  .4194,  with  the  result  that  the 
direct  path  variances  were  then  less  than  those  for  the 
folded  path,  as  expected. 

A  plot  of  the  average  ratios 
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for  this  experiment  is  presented  in  Figure  16.  The  variance 
ratio  for  each  displacement  was  obtained  by  first  determining 
the  ratio  of  each  variance  pair  in  the  run,  then  averaging 


these  together.  Very  few  variance  pairs  were  excluded,  only 
those  having  variance  ratio  values  an  order  of  magnitude 
different  than  the  variance  ratio  mean,  and  those  whose 
variances  were  inverted.  The  solid  line  represents  the 
theoretical  relationship  advanced  by  Ze'evi  for  the  case  of 
Lq  =  20  cm,  which  seemed  to  provide  the  best  fit  of  the 
experimental  data  to  the  theoretical  curve.  We  were  relatively 
free  to  scale  the  outer  turbulence  LQ  as  desired,  as  long  as 
it  did  not  exceed  the  height  of  the  optics  above  the  ground, 
which  was  90  cm. 

Although  the  average  ratios  plotted  in  Figure  16  do  not 
conform  closely  to  the  theoretical  prediction,  several 
reasons  for  this  disparity  are  possible.  First,  the  technique 
of  projecting  a  plane  wave  front  for  the  direct  path  source 
and  then  converting  the  data  to  compare  with  the  folded  path 
intensity  variance  may  not  have  been  strictly  valid.  More¬ 
over,  the  collimating  lens  introduced  a  pronounced  Fresnel 
zone  structure  in  the  direct  path  spot,  which  could  have 
introduced  error  into  the  variance  values.  Second,  the 
overall  strength  of  the  scintillation  may  have  been  too 
low  on  that  particular  evening  for  consistent  data  to  be 
recorded,  given  the  cumulative  effect  of  the  system  errors 
inherent  in  our  equipment.  It  can  be  seen  from  the  data 
that  the  average  strength  of  the  scintillation  that  night 
was  within  the  range  for  weak  turbulence,  but  remained  at 
the  low  end  throughout  the  night. 
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D.  SECOND  OUTSIDE  EXPERIMENT 

In  order  to  correct  for  the  possible  deficiencies  in  our 
first  experiment,  we  introduced  two  changes  to  both  improve 
the  system  and  increase  the  strength  of  the  scintillation. 

We  removed  the  collimating  lens  pair  to  eliminate  the  need 
to  convert  the  direct  path  variance  to  a  spherical  wave 
scintillation  variance,  and  allowed  the  direct  path  laser 
pulses  to  expand  naturally  upon  leaving  the  laser  and  passing 
through  the  splitter.  Then  we  moved  the  entire  experiment 
to  a  new  location  on  the  roof  of  Spanagel  Hall.  A  40  meter 
range  was  chosen  that  was  free  of  any  ventillation  ducts  or 
other  structures  that  might  change  the  uniform  characteristics 
of  the  path.  We  expected  that  uniform  heating  of  the  red 
bricks  on  the  roof  during  the  day  would  result  in  greater 
scintillation  values  at  night  than  at  our  original  location. 

During  the  second  experiment,  data  was  sampled  for  13 
different  values  of  D  from  D  =  10  to  D  =  110  cm,  in  intervals 
of  10  cm.  Two  additional  displacements  of  D  =  5,  and  D  =  150  cm 
were  also  used.  Once  again,  the  data  was  taken  at  night, 
and  no  problems  were  experienced  in  system  alignment  or 
operation.  The  data  was  collected  over  a  period  of  eight 
hours  on  13  November  between  8:30  PM  and  4:50  AM,  and  the 
weather  was  clear,  cool,  and  with  a  slight  wind  present. 

The  increase  in  scintillation  strength  on  the  roof  is 
definitely  present  in  the  increased  values  of  both  the  direct 
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path  and  folded  path  variance.  In  addition,  the  removal 
of  the  collimating  lens  pair  eliminated  the  need  for 
converting  of  the  direct  path  data.  In  general,  the  range 
of  average  ratios  of  the  variances  is  much  greater  too,  and 
a  plot  of  the  ratio 


JVTF) 

V*Cl7 


versus  y0  (  f-  ) 


for  all  displacements  is  presented  in  Figure  17.  As  in  the 
results  from  the  first  experiment,  the  solid  line  represents 
the  proposed  relationship  advanced  by  Ze'evi.  An  outer 
turbulence  scale  of  LQ  *  20  cm  was  chosen  as  the  best  fit  of 
the  data  to  the  theoretical  prediction.  Once  again,  however, 
the  correlation  was  not  as  close  as  had  been  hoped. 

The  points  plotted  in  Figure  17  are  labeled  for  the 
order  in  which  they  were  taken  during  the  night.  Upon 
closer  examination,  it  appeared  that  perhaps  some  relation¬ 
ship  existed  between  when  the  data  for  a  particular  displacement 
was  collected  and  the  resulting  value  of  the  variance  ratio. 

This  lead  us  to  the  conclusion  that  perhaps  the  outer 
turbulence  scale  of  the  range  was  changing  during  time  required 
to  sample  the  scintillation  at  all  values  of  displacement. 

This  would  make  it  impossible  to  reconcile  all  our  variance 
ratios  to  one  plot  with  a  single  characteristic  value  of  Lq. 
Therefore,  we  decided  to  try  to  reduce  the  total  time  required 
to  sample  the  scintillation  at  a  significant  number  of  values 
of  the  displacement  D,  in  order  to  try  and  construct  the 


Figure  17.  F.xper 


first  part  of  Ze'evi's  curve  experimentally.  It  was  also 
possible  that  even  this  would  not  help  because  the  turbulence 
may  not  be  uniform  throughout  the  length  of  the  range. 

E.  THIRD  OUTSIDE  EXPERIMENT 

In  this  our  third  experiment,  we  set  out  to  collect  as 
many  variance  ratios  as  possible  in  the  shortest  time, 
beginning  with  small  values  of  D  and  working  outward  from 
the  detector.  In  this  way  we  hoped  to  correlate  our  data 
with  the  most  critical  part  of  Ze'evi’s  theoretical  curve. 

To  do  this  required  several  changes  in  our  scintillation 
measuring  technique,  the  most  important  of  which  was  to 
alter  the  computer  program  in  the  HP-9825A  to  record  fewer 
data  pairs  and  store  them  automatically  on  magnetic  tape 
while  aligning  the  system  for  a  new  value  of  displacement. 

We  chose  to  record  4  data  pairs  at  each  value  of  D,  and 
we  arbitrarily  set  a  time  limit  of  30  minutes  to  collect 
as  much  data  as  possible.  Additionally,  we  made  a  minor 
change  to  the  folded  path  source  configuration,  arranging 
the  laser,  Pockels  cell,  and  polarizer  on  a  single  base 
that  could  be  slid  as  a  unit  away  from  the  detector, 
simplifying  and  shortening  realignment  time  at  each  successive 
value  of  D.  Previously  the  components  for  the  folded  path 
source  had  been  independently  mounted  and  aligned. 

During  this  experiment,  we  were  not  able  to  collect 
the  full  range  of  13  points  as  we  had  done  earlier,  but 


were  able  to  collect  five  variance  ratios  in  the  first  28 
minutes.  Setup  and  alignment  for  each  value  of  D  was  very 
rapid,  requiring  about  4  minutes  at  each  position,  and  the 
alignment  of  each  spot  from  both  sources  was  checked  visually 
each  time  to  ensure  that  both  remained  centered  on  the 
detector  aperture.  The  data  was  recorded  during  the  night 
of  21  November,  beginning  two  hours  after  sunset.  The 
actual  recording  of  the  4  scintillation  pairs  on  in  the 
NIC-80  required  less  than  a  minute,  compared  to  4  minutes 
in  our  two  previous  experiments  when  we  were  collecting 
16  scintillation  pairs  at  each  value  of  D. 

More  scintillation  was  experienced  during  this  experi¬ 
ment  than  on  the  other  two  nights.  A  plot  of  the  ratio 

-JliiiZ!  versus  y„  (  ®  ) 

VCL)  0  Lo 

for  the  five  displacements  sampled  is  shown  in  Figure  18. 

An  outer  turbulence  scale  of  L  =20  cm  was  chosen  as  a 

o 

best  fit  to  the  theoretical  prediction;  however,  the 
correlation  was  again  poor. 


Figure  18. 


IV.  CONCLUSIONS 


A.  FINDINGS 

The  findings  of  this  work,  based  upon  our  experimental 
results  and  analysis  of  the  problem,  are  as  follows: 

1.  The  system  as  employed  in  this  experiment  is  suitabl 
for  the  comparison  of  directed  and  folded  path  scintillation 
in  the  weak  turbulence  regime. 

2.  Based  upon  the  results  of  our  chi-square  analysis 
of  the  log  normal  data,  both  the  direct  and  folded  path 
intensities  are  log  normally  distributed  above  o^2  =  0.9. 
Below  this  value,  both  distributions  became  severely 
weighted  about  the  mean  and  departed  markedly  from  a 
Gaussian  distribution. 

3.  The  scattering  of  the  variance  ratios  that  we 
measured  at  various  values  of  displacement  does  not  support 
the  theoretical  work  by  Ze'evi. 

B .  SUMMARY 

As  we  worked  with  the  equipment  and  came  to  understand 
the  problem  better,  we  realized  that  the  business  of 
accurately  measuring  direct  and  folded  path  scintillation 
was  rich  in  potentially  puzzling  effects,  hard  to  identify 
and  even  harder  to  eliminate.  Even  if  all  other  effects 
are  understood,  the  very  nature  of  performing  an  experiment 


outside  the  laboratory  where  you  are  subject  to  wind- induced 
vibrations  of  your  optics  or  refractive  inhomogeneities 
in  the  atmosphere,  limits  your  ability  to  obtain  precise 
data.  When  we  began  this  experiment,  we  set  for  ourselves 
the  following  three  goals: 

1.  Construct  a  system  that  would  allow  us  to  simultaneously 
record  the  scintillation  over  both  a  direct  and  a  folded 
path. 

2.  After  recording  the  scintillation,  calculate  the  scintil¬ 
lation  strength  for  each  path  and  obtain  a  measure  of 
its  validity  based  upon  the  accuracy  of  a  log  normal 
assumption  distribution. 

3.  Compare  our  experimental  results  of  the  ratios  of  folded 
path  variance  to  direct  path  variance  at  various  values 
of  displacement  in  an  attempt  to  verify  the  theoretical 
work  of  Dr.  Ze'evi. 

Of  these  three,  the  first  two  were  achieved  and  the  third 
requires  further  work.  We  were  able  to  construct  and  test 
a  system  that  accurately  and  simultaneously  measured  direct 
path  and  folded  path  scintillation.  This  system  was 
successfully  tested  and  produced  data  that  was  in  agreement 
with  the  accepted  log  normal  distribution  for  scintillation. 

As  expected,  folded  path  variance  was  consistently  greater 
than  the  direct  path  variance  during  all  three  outside 
experiments  utilizing  two  different  optical  configurations, 
and  the  range  of  the  variance  ratios  was  between  2  and  7 
in  magnitude.  Applicability  of  the  weak  turbulence 
scintillation  theory  was  valid,  based  upon  the  fact  that 
all  log  intensity  variances  were  well  below  1.2.  Additionally, 


the  statistics  program  that  we  wrote  for  the  HP-9825A 
worked  well,  and  the  chi-square  test  proved  to  be  an 
accurate  measure  of  the  data's  correlation  to  a  Gaussian. 

What  we  were  not  able  to  do  in  our  experiment  was  to 
conclusively  prove  or  disprove  the  folded  path  theory.  As 
can  be  seen  from  the  summary  plots  of  experiments  one  and 
two  depicted  in  Figures  16  and  17,  there  is  a  fairly  wide 
scattering  of  the  variance  ratios  at  all  displacements, 
and  the  entire  collection  of  variance  ratios  from  experi¬ 
ment  three  depicted  in  Figure  18  lies  well  below  the 
theoretical  curve.  Taken  at  face  value,  the  data  does  not 
support  the  theory;  however,  the  wide  range  of  variables  in 
this  experiment  that  must  be  controlled  to  obtain  valid 
results  certainly  leaves  open  the  possibility  that  the 
theory  is  correct  and  that  the  experimental  method  was  not 
precise  enough. 

C.  EXPLANATIONS 

To  answer  the  question  of  what  factors  may  have  been 
responsible  for  the  scattered  nature  of  the  data,  it  is 
necessary  to  look  at  the  data  collection  techniques  that 
we  employed.  The  first  factor  that  may  have  played  a  part 
is  the  quality  of  the  sources  that  we  used.  In  Chapter  II 
we  stated  that  the  correlation  and  variance  of  the  HeNe 
sources  had  been  examined  by  Hodgini  and  found  satisfactory. 
It  is  possible  though  that  with  time  and  usage,  especially 
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considering  the  amount  of  moving  about  from  place  to  place 
that  we  subjected  the  lasers  to,  that  their  internal  stability 
may  have  been  affected.  The  light  from  these  lasers  is 
randomly  polarized,  with  a  source  intensity  variance  of 
approximately  10’4.  If  this  could  have  increased  with  time, 
then  it  is  certainly  possible  that  the  sources  would  have 
contributed  to  the  overall  scintillation  that  we  measured. 

The  second  factor  that  may  have  played  a  part  is  that 
of  changing  atmospheric  conditions,  specifically  the  outer 
turbulence  scale  length  LQ.  This  could  occur  during  the 
collection  of  a  pair  of  path  distributions,  during  the 
collection  of  an  entire  set  of  16  pairs  at  a  given  displace¬ 
ment,  or  during  the  time  required  for  a  set  of  displacement 
values  to  be  measured.  If  any  inhomogeneities  in  the 
turbulence  were  present,  or  if  the  turbulence  was  nonuniform 
throughout  the  length  of  the  folded  path,  then  this  would 
certainly  have  affected  the  data  as  well.  Looking  at  the 
first  two  experiments  together  since  they  both  employed 
the  same  general  technique  of  long  periods  devoted  to  the 
collection  of  the  data,  it  seems  reasonable  to  assume  that 
non-steady  state  atmospheric  conditions  were  probably  the 
dominant  effect  involved  in  destabilizing  the  data.  This 
is  partially  borne  out  by  the  results  of  our  third  experi¬ 
ment,  in  which  the  variance  ratios  resulting  from  a  much 
shorter  period  of  time  are  significantly  more  consistent. 
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Taken  in  the  aggregate,  the  data  is  extremely  inconclusive. 
(  No  general  factor  or  pattern  of  the  variance  ratios  emerged. 

One  fact  does  stand  out  however,  and  that  is  this:  future 
work  in  this  area  will  absolutely  have  to  control  the  uni¬ 
formity  and  magnitude  of  the  turbulence  in  the  path  geometries 
to  accurately  verify  the  folded  path  theory.  Ideally,  all 
12  or  13  displacements  of  the  folded  path  should  be  measured 
simultaneously  in  a  wind  tunnel  of  controlled,  uniform 
turbulence  in  order  to  be  absolutely  certain  that  the  same 
factors  present  in  one  measurement  will  be  present  in  all. 

Any  past  attempts  to  correlate  folded  path  scintillation 
with  that  of  direct  path  scintillation,  and  which  did  not 
take  into  account  this  variable  nature  of  the  turbulence 
must  be  considered  suspect. 

D.  FUTURE  WORK 

This  experiment  needs  to  be  repeated  using  better  sources 
and  an  optical  configuration  that  will  allow  the  nearly 
simultaneous  measurement  of  the  folded  path  scintillation 
for  several  values  of  displacement  under  uniform  conditions. 

If  better  and  higher  power  sources  were  employed,  the  source 
stability  problem  could  be  eliminated,  and  the  experiment 
could  be  performed  during  the  daytime  utilizing  a  narrow 
bandpass  filter  in  front  of  the  receiver  and  neutral  density 
filters  to  control  the  intensity  of  the  direct  path  signal. 
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A  broader  data  base  should  be  constructed,  including 
comparison  of  direct  and  folded  path  scintillation  at  various 
heights  above  the  ground.  Further  investigation  into  the 
effects  of  aperture  averaging,  different  weighting  functions, 
and  target  interaction  as  they  relate  to  direct  and  folded 
path  scintillation  is  needed.  The  effect  on  direct  versus 
folded  path  correlation  if  real  world  variables  such  as 
vibration  are  included  needs  to  be  investigated  as  well. 

These  areas  must  all  be  addressed  if  the  theory  is  to  be 
usable  to  predict  optical  system  performance  under  realistic 
field  conditons. 


APPENDIX  A 

NIC- 80  COMPUTER  PROGRAM 

On  the  pages  that  follow  is  listed  the  computer  program 
we  used  for  the  Nicolet  Instrument  Corporation  Model  NIC-80. 
The  program  is  written  in  the  assembly  code  for  the  machine 
which  is  given  in  the  operating  manual,  and  is  compiled  on 
and  loaded  from  the  HP-9825A.  The  9825A  program  to  perform 
these  functions  is  an  in-house  program  written  by  Professor 
Edmund  Milne  at  the  Naval  Postgraduate  School. 
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257 

1000240 

JMP @  240 

260 

330  3  77 

M-A 

A  (  377 

26  1 

5006 

LAS 

6 

2o  2 

214001 

TAC 

26  3 

3135670 

MPO 

AM@  16  70 

Ji 

264 

10  5000 

VDS 

0 

26  5 

4G13 

T¥S 

266 

1000240 

JMP  («  240 

267 

350 

JMP 

350 

2  7u 

2111700 

MEM 

A  1700 

SO  0 

271 

2  40  5o  77 

ACC 

a  16  7  7 

I 

272 

110400 

MEM 

A  (  400 

273 

240567b 

ACC 

M  16  76 

J 

274 

3111677 

MEM 

Ais?  loll 

I 

275 

2000015 

JMS 

15 

276 

2125677 

MPO 

M  16  7  7 

I 

277 

2707676 

MMO 

MZ  lo  76 

J 

77 


Scintillation  File  #  9  Tracn  #  0 


Address 

Content 

Mnemonic 

300 

274 

JMP 

274 

301 

1000267 

JMP@ 

267 

302 

3001723 

JMS  (3 

1723  RES 

30  3 

2000240 

JMS 

240 

304 

6464 

HSR 

30  5 

30  3 

JMP 

303 

30  6 

44463 

RHS 

307 

405160 

EZA 

310 

0 

JMP 

0 

311 

110003 

MEM 

A(  3 

312 

3405727 

ACC 

M(e  1727 

DSK 

313 

2111700 

MEM 

A  1700 

BOO 

314 

3405726 

ACC 

MIS  1726 

DSA 

315 

3001725 

JMS  is 

!  1725  DSP 

316 

44463 

RHS 

317 

405160 

EZA 

320 

0 

JMP 

0 

321 

472001 

A-M 

AZ  (  1 

322 

2162000 

ZER 

Z  0 

32  3 

302 

JMP 

302 

324 

472001 

A-M 

AZ  (  1 

325 

2162000 

ZER 

Z  0 

326 

347 

JMP 

347 

327 

462001 

A-M 

Z(  1 

330 

350 

JMP 

350 

331 

2000267 

JMS 

267 

3  32 

111000 

MEM 

A (  luOO 

333 

2405677 

ACC 

M  16  77 

I 

334 

110400 

MEM 

A  (  400 

335 

2405676 

ACC 

M  16  76 

0 

336 

64  64 

HSR 

3  37 

336 

JMP 

336 

340 

4446  3 

RHS 

341 

5004 

LAS 

4 

342 

3405677 

ACC 

M  16  7  7 

I 

343 

2  125677 

MPO 

M  lo  7  7 

X 

344 

2  7  0  7  o  76 

MMO 

HZ  1676  j 

345 

3  36 

vj  MP 

33o 

346 

0 

omP 

0 

347 

2u0026  7 

JMS 

267 

350 

2 1 1 1 6  76 

MEM 

A  lo  76 

J 

351 

4  7 luOO 

A-M 

A(  lo  U  o 

352 

3405726 

ACC 

M  i?  1726 

CS  A 

353 

3u  o  1 7  2  5 

J  MS 

4  1 7  c.  5  u 

SP 

354 

4446  3 

RHS 

355 

12u 

uMP 

12u 

78 


APPENDIX  B 


HEWLETT-PACKARD  982 5 A  CALCULATOR  PROGRAMS 

Two  programs  for  the  HP-9825A  are  listed  in  this 
appendix.  The  first  program  does  not  give  the  operator  the 
option  of  storing  the  data  run  Cat  a  specific  displacement) 
on  magnetic  tape  before  advancing  to  the  next  displacement. 
The  second  program  does  contain  this  option,  as  well  as  the 
choice  of  how  many  data  pairs  (pairs  of  direct  and  folded 
path  scintillation  distributions)  to  record  at  each 
displacement.  Note,  however,  that  the  number  of  data  pairs 
selected  in  this  program  MUST  correspond  to  the  number  of 
data  pairs  specified  in  steps  110  and  221  of  the  NIC-80 


program. 


FIRST  PROGRAM 


0 

1 

2 

3 

4 

5 

6 

7 

8 
9 


■"SPEER/PARKER  SCINTILLATION  PROGRAM  8/8  2": 
asp  "SCINT  PROGRAM:  push  contin";stp 
dim  D&t  13]  rMS (2  ,6  0]  ,0(10] 

aim  A$(3S,4G]  ,BS(7]  ,P.[256]  ,G[256]  ;716-F;701-E 

aim  S{2  56J  ;aev  "E"  ,E,  "F"  ,F;on  err  “erl" 

dim  N$(12,3] ,R[256] , I [ 2 56 ] 

fxa  0;iaf  2,A$;mdec?"  meters"*B$[l,7] 

time  250? on  err  "erl" 

it  ras ( 701)  #17 ;asp  "Printer  is  not  ready";stp 
6.  228e-7-X;33.169-r0;l-Q(l]  ;  5-C(2]  ?  S-Q [ 3)  ;13-G(4]  ?17-Q(5] 


10 

J1 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 
27 
26 

29 

30 

31 

32 

33 
24 
35 
3b 
37 
3e 

39 

40 

41 

42 
4  3 
h4 
45 
*»6 


wtb  E, 27, 69, 27, 84, 27, 76, 13, 32, 32, 32, 32, 32, 32, 32; 21-0(6] 
wtb  E, 32, 32, 32, 32, 32  , 27, 77 

25-0(7]  ?  2S-0(6J  ;.33-0 [  9]  ?  37-0(10]  ;.if  ras(2)  #296?cll  'er2  ' 
0-L ;  l/r0-r  25;tor  J  =  1  to  256;  (J+L)  r25/2-S. [J ]  ;J-L; next  J 
"JAN" -NS  ( 1]  ;  "  F  EB  "  -N  $  (  2  ]  ;  "MAR"*N$[  3] 

"APR"  -N  §  (4  ]  ;  "MAY"  *N$  (  5  ]  ;"JUN"*N$(  6]  ?  "  JUL"-N$  (  7] 

"  A0G  "  -N  $  ( 8  ]  ?  "  SEP  "  *N  $  [  9] 

"GCT"-N8 ( 10 ] ;"NCv"-N$(ll] ; "DEC"-N$( 12] 
on  err  "er2";wtb  2,l;gto  "scint" 

"er2" : beep ;csp  "NIC  80  DOWN" ;wait  250;asp  "  " ;jmp  0 
"erl" :Deep?asp  "error  “ ,char ( rom) ,ern, "  line  “ ,erl 
wait  1000 

if  ern=66;wrt  "E"  , "Div isicn  by  zero  occurea  in  line", erl 
stp  ;gto  ,20 

"input":ent  "Enter  the  number  of  tne  montn",N 
ent  "Enter  tne  aay"  ,T 

N§[N,1  ,3]-D$(l,3];".  "-08(4,5]  ;  str  (.T)-C$(  6 ,8] 

"  1962"-D8(9  ,12] 

rlt  3;asp  "Enter  waveleng  tn "  , X;  ent  " "  ,x,*txa  0 
ent  "Enter  run  ccmroer ts"  ,M$  ( 1 , 1, 60  ] 
ent  "Enter  weatner  comments" ,M$ [2 ,1 ,60] ;ret 
"scint" :on  err  "erl" ; time  30000 

ent  "Enter  u  tor  cal.  cr  0  rcr  cate", M; it  M#0?gtc  "cal" 
ent  "Enter  range"  ,R 

ell  ’'input'; asp  "set  tcp  ot  torm  on  pr inter "; stp 
wtb  :2,M;wtb  E, 27, 69, 27, 84, 27, 76, 13, 32, 32, 32, 22, 32, 27, 77 
tor  1=1  tc  16 ;asp  I, "Scint  accumulation" 
rab(2)-J;if  J#I;wtb  2, 1;  beep;  asp  "OU.I  CF  SYNC"  ;  stp 
wtb  2, I; next  I 

ent  "Enter  0=acort,  l  =  cutput"  ,B;.ir  B  =  u;wtD  2,£;gtc  0 

tmt  1 , cl3  ,  /  ,  c6 0  ,  /  ,  c6 0  ,/  ,  "Wave  ler.g  tn=  "  ,el0 . 2  ,  c7 , /  ,/ ,/ 

tmt  2,"  t  pts  Sigma  Cm  Square",/ 

rmt  3,r4.0,":",f5.0,x,fl0.4,x,  lelO . 2 

rmt  4  ,r4.0 ," ,t5.C,x,f 10. 4 ,x , lei 0.2, "**" 

rmt  6,"Fixea  range="  ,  to.  0  ,c7  ,/ 


wrt  "£.1"  ,D$,M$(  1,1,60]  ,M$  {  2 , 1 , 6G]  ,X,E$ 

wrt  " E" , "Scintillation  measurements  or  Sigma" 


f 


wrt  "E "  ,"16384  samples  for  each  measurement" 

wrt  "E  .  6  "  ,  R  ,B$  ' 

wrt  "E.2";gto  50 

for  1*1  to  16;wtb  2,I;rdbi2)*B 

if  B#0; beep; dsp  "NOT  IN  SYNC"; wait  2000 ;gto  0 

for  J*1  to  256;  1024*rdb<2) +rdb(2)-R[J]  ;next  J 

wtb  2,l;for  J.*l  to  256;  1024*rdb(  2)  +  rdb(  2) -I  [  J] 

next  J;0*r20*r21 

for  J=1  to  256;r20+R[J] +r  20;  r2  1+ 1  [  J] +r21  ;next  J;0*B*A 

csp  "Statistics  Calculations" 

for  J*1  to  256;A+R[J]  *S.[J] -*A;B+I[J]  *S[J]«-B 

next  J ; A/r 20*A; B/r 21+B ; 0*S*V 

for  J=1  tc  256; S+R [ J] (S(J ] -A) (S [ J ] -A) +S 

V+L[J] (S[J] -B) (£[ J] -B) *V;next  J 

V  ( Vr21)  ♦V;  y  (S/r  20)  ♦S  ;S+rl8;  r2  0+Z;  r25/2+r6 

for  K= 1  tc  2; 0*r 17; f or  J=1  tc  256 ;S[J] -r6*r7 

S(J]+r6-r8;  (r7-A) /r  18*1Q0-t9;  ( r8-A)/rl8*100-*-rl0 

f  rc(abs(r5) )  *r30;  frc(abs(rlO)  )  *r26  ;r9/10*r9;r  10/10+r  10 

for  Y*1  to  2 

int(abs(r  9)  )-»rll;int(abs(rlO)  )*rl2;  if  r  9<Q;  -  r  ll-»-r  11 

if  r 10<0; ^rl2-r 12 

if  rll>36;  3a-rll+r9-r9;.38-rll 

if  r ll<-38 ;  \t 11  +  38)  ( -1) +r9-r 9; -38*r 1 1 

it  rl2>38;3a-rl2+rl0-rl0;  38-rl2 

if  rl2<-38;-(rl2+38) +r 10+rl0;  -38*rl2 

int(frc(abs(rS))*10+l)*rl3;int(frc(  abs  (rlO)  )  *10+1)  -»-r  14 

Ci  [r  13]  +r  13  ;w[  r  14]  *rl4 

if  rll>0;stt  (A$[rll,rl3,rl3+3]  )-*-rl5 

it  rll*0;stf  (A$[39,rl3  ,rl3+3]  )*rl5;if  r  9<0;  l-rl5*rl5 
if  rll<0;abs(rll)-*-rll; 1-stf (A$(rll,rl3,rl3+3])*rl5 
if  rl2>0;str  <A$(r  12  ,rl4  ,rl4  +  3J  )  ♦rl6 

it  r  12=0  ;  stf  <A$  [  39  ,  r  14  , r  14  +  3]  )  ♦rl6  ;  if  r  L0<0 ;  1-r  16 ♦r  16 
if  rl2<0;abs(r]2)-*rl2;  1-stf  (A8[rl2,rl4,rl4+3])-*-rl6 
if  Y*l;r  16  +  r28;r  j5-*r27;ytc  82 
gto  66 

if  r  5<0;  r9- .  l*r  5  ;gtc  84 
r  9+ . I*r9 

it  r 1G<0;  rlO- „l*r 10; gtc  86 
r  10+  .  1  *c  lu 

next  Y;  (cl5-r27)*r20+r27*rl5;(rld-r2B)*r2b+r2tt*ri6 

prnci<Z*(rIo-rl5)  ,U)-P[j] 

it  P  (u  j  >u;  (  R[  J  ]  -P.  (J  ]  )~2/P{GJ+rI7-*-rl7 

next  J ; it  K=2;gtc  92 

it  r 17 >2b 7 . 86 ; wr t  ”  £.  4"  ,1  ,'L  ,S  ,  rl 7  ;cll  ‘  'plot ' ;  g tc  94 
wrt  ”  t.  3"  ,1  ,2  ,S  ,  r  17;  ell  ’ 'Piet '  ;g  tc  94 
if  rl7>287.BO;wr  t  "E.  4"  ,1  ,2  ,v,  rl7  ;c  11  'Plot gtc  95 
wrt  "E.  3",I  ,2,  V,  rl7;cll  'Plot';gtc  9  5 
051 


94:  r21-*-Z  ;B*A;  V*rl8;ara  I+R 

95:  next  K 

96:  wtb  2#l;next  I 

97:  wtb  E,12;gto  7 

96:  "cal"  :wtb  2,1 

99:  ent  "O-abcr t, L-aisplay " ,D; wtb  2,D;if  D=Q;gto  0 

100:  ent  "Q-abor t, 1-reset , 2-c  out, 3-bin", D 

101:  if  D=l;wtb  2,D;jir.p  -2 

102:  if  D=0;wtb  2,D;gto  0 

103:  if  D  =  2;wtb  2,l;gtc  "c  cut" 

304:  if  D  =  3;wtD  2,D;gto  "bin" 

105:  "c  out" :dsp  "DON'T  USE  W/O  DIGITAL  VOL1ME1ER! " 

106:  stp  ;dsp  "Set  VI, push  continue" ; stp 

107:  rea  722,K;red  722,K;asp  "Set  V2,pusn  continue"; stp 

108:  red  722,L;rea  722, L 

109:  tor  1=1  to  256;  1024*rab  (  2) +rab(  2) -»-A 

110:  if  AMO;  txa  3;prt  A+I/1000;rxa  0 

111:  next  I 

112:  ent  "Enter  XI"  ,r 2  ."Enter  X2",r3 
113:  (c3-r2)/ ( In  (L).-ln  (K).)-r0;gto  0 

114:  "bin " :  0-*A ; f or  J  =  1  to  2 56 ;  1024  * rab( 2) + rab ( 2) -I fu ] ; next  J 
115:  ir,  l  n  ( I  (  *  ]  )  +A 

116:  for  0=1  tc  256; int ( 32000*^/1 (0 ] )*E; wto  2,E;next  0 
117:  gee  7 

11b:  "Plct":ent  "PLOT  type  1, continue  0"  ,r22; if  r22#l;gtc  3  4  3 
119:  csiz  1. 7 , 1. 7  ,  .7  1 , 0 

120:  fxa  0  ;scl  .-1 . 3 , 6 . 9  P-15G  ,  1600  ;  pi  t  6.2,1550,1 
121:  lbl  "Plot  of  " ;  lbl  K 
122:  lbl  "  Aay";cplt  -14, -1 

123:  lbl  "Scant  illation  Run  fr " ; lcl  I;plt  6  .2, 1455, 1; iol  DS 

124:  cplt  -13,-l;lcl  MS  [  1 , 1 ,60]  ;  rxa  5 

125:  pit  6.2, 1360, l;lol  ( 2 , 1 , 60] ; pi t  6.2,1313,1 

126:  lbl  "Exper.  Sigir.a  " ;  lcl  rl  6 

127:  pit  6.2,1265,1 

12b:  lbl  "Exper.  Mean  ";lb 1  A;txd  3;plt  6.2,1217,1 

129:  lbl  "Cm  Square  ";lcl  rl7;txc  0;pen 

130:  pit  0,0,1; pit  6. 7, 0,2= pit  0,0,1; pit  0,lb00,2;pen 

131:  ter  x=u  tc  8  cy  2 ;  pi  t  X,-9,l;cplt  -.5,-1.5;lcl  X;r.ext  X 

132:  0  +  Y 

133:  pit.  -  .  1 , Y , 1 ;cpl t  -2,-.3;lci  Y  ;  Y+5Gu-»Y  ;  ir  Y=2000  ;gtc  135 
134:  g  tc  133 

135:  pit  3,-15C,l;lcl  "Log  I  r.tensi  ty  "  ;  pi  t  -  1,500  ,1 
136:  csiz  1 . 7  , 1. 7 , . 7 1 , 90; lbl  "Nc.  or  Counts" 

137:  csiz  1  .  7  ,  1 .  7  ,  . 7  1 , 0;  ter  J=1  tc  256 
13c:  pit  £  [  J  ]  ,  R  lo  ]  ^-2;  per.;  nex  t  J;pen;pit  0,0, -2 
139:  tor  J=1  tc  256;  pit  £  [0  ]  ,  P  ( j  ]  ;  r.ex  t  J 
140:  pit  o.2, 1169,  l;lci  ".  Exper.  Curve" 

141:  pit  b. 2, 1122,1;  lbl  Cr.i  £u.  Preuicticn" 

1*42  :  pit  6  ,  1500  ,  i 


SECOND  PROGRAM 


0:  "SPEER/PARKER  SCINTILLATION  PROGRAM  8/82": 

1:  asp  "SC INI  PROGRAM:  pusn  contin";stp 
2:  aim  DS(13]  ,MS(2, 601  rX  ,BS(7]  ,R  ,0(10] 

3:  aim  A$ ( 39 ,40]  ,P [ 256] ,G { 256] ; 716*F; 701*E 
4:  dim  S(256J ;aev  "E"  ,£, MF" ,F;on  err  "erl" 

5:  aim  N$(  12, 3]  ,R[256]  ,1  [256] 

6:  txa  0;lar  11  ,A$;  maec;  "  meters“*B$ ( 1 , 7] 

7:  time  250; on  err  "erl" 

8:  ir  ras ( 70 1) #17; asp  "Printer  is  not  reaay";stp 

9:  6. 328e-7*X; 33. 169*r0; 1*0(1]; 5*012] ; 9*0(3] ;13*C(4] ; 17*0(5] 

10:  wtb  E, 27, 69, 27, 84, 27, 76,12, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32 

11:  wtD  E, 32, 32, 27, 77 

12:  21*0(6]  ,*:25*G(7]  ;  2S*0(8] 

13:  33*0(9] ;37*O(10]; it  ras < 2) #296 ;cll  'er2  ' 

14:  U*L  ;  l/rO*r  25;tor  J  =  1  to  256;  (J+L)  r25/2*S[J]  ;J*L;next  J 
15:  "JAN  "*N$  [  1]  ;  "  F  EB  "  *N  8  (  2  ]  ;  "MAR"*N$[ 3] 

16:  "APR"  *N$  [ 4 ]  ;"KAY"*N$(5]  ;  "  JUN" *N$  (  6]  ;  " JUL"*N$. [ 7] 

17:  "AUG" *N$ [ 8 ]  ; " SEP " *N $  (  9] 

18:  "CCT"*N$(  10  ]  ;"NCV"*N$[  11]  ;  "DEC  "*N$.  [  12  ] 

19:  on  err  Mer2";wtb  2,l;gto  "scint" 

20:  "er2'":beep;asp  "NIC  80  DOWN";wait  250;asp  "  "jjmp  0 
21:  "erl" :beep;asp  "error  "  ,cbar ( rom)  , ern, "  line  "  ,erl 
22:  wait  1000 

23:  it  ern=66;wrt  "E"  ."Division  by  zero  in  line" ,er l;stp 
24:  gtc  21 

25:  "input"  :ent  "Enter  tne  number  or  tne  montn",N 
26:  ent  "Enter  the  aay"  ,T 

27:  NS [N , 1 , 3  J  *DS ( 1  , 3  J ; " .  "  *D6  [4,5]  ;str('L)*D8[6,6] 

28:  "  1982"*D$(9  ,13] 

29:  tit  3;asp  "Enter  wavelengtn"  ,X;  ent  "",X;txa  C 
30:  ent  "Enter  run  comments"  ,M$ [ 1 , 1 , 60 ] 

31:  ent  "Enter  weatner  comments"  ,M$ ( 2 ,1 , 60] 

32:  ent  "Enter  pair  counter  ",  r3b  ;Q*r31 
33:  ent  "Enter  record  tile" ,r32; r32*r34 
34:  ent  "Enter  tracK  num.ber  "  , r35 ;  ret 
35:  "scint ":on  err  "erl"; time  306UO 

36:  ert  "Enter  1  tor  cal.  or  0  tor  aata",iy;ir  M#0;gto  "cal" 
37:  ent  "Enter  range"  ,R 

38:  c1  ';nput';csp  "set  top  ct  term  cn  printer" ;stp 

39;  w ,M ;  w to  E  , 2 7 , 6 9  ,2 7  ,84  ,27 , 76  ,12 , 32 , 32 , 32 , 32 , 32 , 32 , 32 

40:  wtb  £,32,32,32,32,27,77 

41:  tor  1=1  tc  r36;asp  I, "Scmt  accumulation" 

42:  rab(2)*J;it  J#I;wtb  2,I;beep;asp  "CUT  Cf  SYNC";Stp 
43:  wtb  2,1; r31+l*r 31; nex t  I 

44:  esp  "Recording  Cata";trK  r35;ctg  ; ret  r 32 ,D $ ,M$ , X , E$  ,R 
45:  r32+l*r32 

46:  tor  1=1  tc  r36;wtb  2  ,1 ;  rob  (  2)  *B 


47:  if  B#0 ;beep;asp  "NOT  IN  SYNC";wait  2000;gto  0 

48:  for  J=1  to  256 ; 1024* rdb{ 2) +rcb(2) ♦RfJ] ;next  J 

49:  wtb  2,l;for  J=1  to  256 ; 1024*rab( 2) +  rab( 2) ♦L [ J ] ; next  J 

50:  ref  r32,Rt*.] ,1 [*]; r22+l*r 32; wtb  2,l;next  I 

51:  ent  "Enter  0  next  run,l  for  calc",r33 

52:  if  r33=0;dsp  "READY:  Push  Cont.";stp  ;gto  38 

53:  for  M=1  to  r31;ent  "Enter  start  file",r34 

54:  ent  "Enter  track  number ", r35; trx  r35 

55:  laf  r 34 ,D$ ,M$ ,X,B$ ,R; r 34+l^r 34 

56:  fmt  1 , cl3 ,/,c6 0 ,/ , c6 0 ,/ , "fta veleng th=  "  ,el0 . 2  , c7  ,/ ,/ ,/ 

57:  fmt  2,"  #  pts  Sigma  Chi  Square",/ 

58:  fmt  3  ,  f  4 . 0  ,  "  :  "  ,  f  5 . 0  ,x  ,f  10 . 4  ,x  ,  lelO .  2 
59:  fmt  4  ,  t4. 0  ,  "  :  "  ,t 5 . 0  , x  ,t  10 . 4  ,x  ,  lelO .  2  , "  **" 

60:  fmt  6,"Fixea  range®"  ,f  6. 0  ,c7  ,/ 

61:  wrt  "E.  1"  ,D$,M$[1, 1,60]  ,K$[2  ,1,60]  ,X,B$ 

62:  wrt  "E " , "Sc intillation  measurements  of  Sigma” 

63:  wrt  "E", "16384  samples  for  each  measurement" 

64:  wrt  "E .  6"  ,  R  ,  fi$ 

65:  wrt  "E  .2" 

66:  tor  1=1  to  r36;laf  r  3  4  ,  R  [  *  ]  ,  I  [  *  ] 

67:  0*r20*r21 

68:  tor  J=1  tc  256 ;  r20+R  [  J  ]  ♦r  20  ;  r2 1+ 1  [  J]  ♦r  21 ;  next  J;0^E-*-A 

69:  osp  "Statistics  Calculations" ; for  J=1  to  256 ; A+R [J ] *S [ J ] ♦A 

70:  B+L [ J] *S[J] ♦£; next  J ;A/r20-A ;B/r21-B;0*S-V 

71:  tor  J=1  to  256 ;S+R [ J ] (S { J ] -A) ( S [ J ] -A) ♦S 

72:  V  +  L[J]  (&[J]-E)  (S(J].-B)  ♦Vjnext  J 

7  3:  V  <  v: / r  2 1 )  ♦\/;y  (S/r2C)  ♦S;S-rl&;  r2 0*Z ; r25/2*r6 
74:  tor  K= 1  tc  2;0*rl7;tor  0  =  1  to  256; & ( J ] -r6 *r 7 
75:  &[J]+r6-r8; ( r7-A) /rl8*100-r9; (r 8-A)/rl8*100-rl0 

76:  t  rc  (abs  (r  i) )  ♦r30;  tr  c(abs  (rl  0)  )  *r26  ;  rS/10*r9  ;r  10/ 10 *r 10 
77:  tor  Y=1  tc  2 

78:  int (abs (r 9) ) ♦  rll; int (abs(rlO) ) ♦rl2;  if  r S<0; -rll+r II 

79:  if  rl0<0; ->rl2^rl2 

80:  if  rll>36  ;  3&-rll+r9*r9;.28*rll 

81:  if  rll<-38;  <r  Jl+38)  {-L)+r9*r  S;.-38  +  rll 

82:  if  rl2>38; 38-r 12+r 10-rl0; 38-rl2 

83:  if  rl2<-36?- »rl2+38)  +r  10^rl0; -26-*>rl2 

84:  int (f  rc (aos (r  £)) *10+1)  +  r  12;  int  ( f  r  c  ( abs  (rlO)  )*10+l)+rl4 

85:  0(rl3]*r 13 ;Q [r 14] *rl4 

8b:  it  rl  1>Q  ;stf  (A$.  [r  11 ,  rl3  ,  rl3+ 3]  )  ♦rlS 

87:  if  rll=0  ;stf  (AS  (  3S  ,  rl3  ,rl  3+ 2]  )  ♦rl  5;  if  r 9<0;  1-r 15^rl5 
88:  it  rll<0; abs(rll) ♦rll; 1-stf (A$[ rl 1 , r 1 3 , r 13+3 ] ) ♦rlS 

8  9:  it  rl2>0;stf  (AS1  r!2  ,rl4  ,r!4  +  3]  )  ♦rlfi 

90:  if  rl2=0;stf  (A$.[39,rl4  ,rl4  +  3]  )  ♦rl6;.it  r 10<0; l-rl6*r 16 
91:  if  r J2<0 ;aos ( rll) ♦rl2; 1-stt ( A$[r 12 , rl4 , rl4+3] ) *c 16 
92:  if  Y=l; r 16*r28; rl5*r 27 ;gtc  94 
93:  gtc  98 
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94:  if  r9<0; rS-. I*r9;gtc  96 
95 :  r  9+.  l-»-r9 

96:  if  rl0<0;rlG-.l+rl0;gtc  98 
97:  rlQ+.l+rlO 

96:  next  i ; (r ]5-r27) * r30+r 27*r 15; ( r!6- c28) *r 26+r28-r 16 
99:  prnci  <Z*  (r  16-r  15)  ,  G)  *P.  (J  ] 

10G:  it  P.(J]  >0;  (R[J]  -P  [J  J  )~2/P.{J]  +rl7*rl7 
101:  next  J;if  K=2;gtc  104 

102:  if  rl7>287.88; wr t  "E  .4"  ,1  ,Z  ,S  , rl7  ;cll  ' 'Plot ' ;  g to  106 

103:  wr t  "E.  3"  .1  ,Z,S , r  17 ;cll  ‘ 'Plot '  ;g  to  I0o 

104:  if  rl7>267.88;wrt  "E  .4"  ,1  ,Z  ,  V  ,rl  7;  c  11  ' 'Plot gtc  107 

105:  wr  t  "E.  3"  ,1  ,Z,V  ,r  17,-cll  'Plot  ";g  to  107 

106:  r2 1+Z ;E+A; V-*>r 18;ara  I+R 

107:  next  K 

106:  r 34+1+r 3 4 ; nex t  I 

109:  wtb  E, 12; next  M 

110:  wtb  E,12;gtc  0 

111:  '"cal"  : wtb  2,1 

112:  ent  "0-abort, 1-oisplay" ,D;wtb  2,D;ir  D*0;gto  0 

113:  ent  "0-abort, 1-reset ,2-c  out, 3-bin", D 

114:  ir  D  =  l;wtb  2  ,D  ;  jmp  .-2 

115:  if  D=0;wtb  2,D;gto  C 

116:  if  D  =  2;wtb  2,l;gto  "c  out" 

117:  it  D  =  3;wtb  2,D;gto  "bin” 

lib:  "C  out"  :asp  "DON'T  USE  ROUTINE  W/O  DIGITAL  VOLTMETER!" 

119:  stp  ;asp  "Set  vl  ,pusr.  continue" ; stp 

12U:  rea  722,K;rec  722,K;asp  "Set  V2,pust  cent inue "; stp 

121:  rea  722,L;rea  722, L 

122:  for  1*1  tc  256; lC24*rab( 2) +rab( 2) ♦A 

123:  it  A>10;txa  3;prt  A+I/lC0O;rxd  0 

124:  next  I 

125:  ent  "Enter  XI"  ,r 2 , " Enter  X2" ,r3 
126:  (  r  3-  r2  )  /  ( In  (L).-ln  ( K).)  ♦rC  ;g  to  0 

127:  "bin" :  U-*A;ror  U=1  to  2  56  ;  1024  *  rab  (  2)  +  rab  (  2)  *1  [  J  ]  ;  next  J 
12  8:  mi  n  ( I  [  *  ]  )  *  A 

129:  ter  0=1  tc  256; int (22uOO*A/I [ J ] ) ♦£; wtb  2, b; next  J 
120:  gto  7 

121:  "Plot":ent  "PLOT  input  1, continue  0"  ,r22; ir  r22*l;gto  j5<* 
132:  csiz  1. 7  , 1 . 7  ,  . 7 1 , 0 

133:  rxc  C  ;scl  .-1 . 3 ,  tt  .  9  ,.-150 , 1600  ;  pi  t  6.2,1550,1 
124:  ltl  "Plot  or  " ;  1  b  1  K;lbl  "  way";.cplt  -14,-1 
125:  ltl  "Scintillation  Run  ;lol  I;plt  6 . 2 , 1455 , 1 ; ltl  D$ 

126:  cplt  -13,-l;lol  [  1  , 1  ,6Cj  ;  r  xa  5 


137: 

pit 

6. 2, 1360,1; lei 

M  $  [  2 

,1,60] ; pit  6.2,1313,1 

128: 

ltl 

"txper.  Sigma 

”;  ltl 

rl  8;  pit  6.2 , 1265  ,  1 

139: 

ltl 

"Exper.  Mean 

"  ;  1 1 1 

A;rxa  3; pit  6.2,1217 

140  : 

ltl 

"Cm  Square 

"  ;  ltl 

a  !  7  ;rxa  C;pen 

141:  pit  0,C,l;plt  8. 7, 0,2;  pit  0,G,l;plt  0,1600, 2;pen 

142:  tor  X»0  to  8  by  2;  pi  t  X  ,-9 ,1  ;cpl  t  -  .5,-1. 5;lbl  X;  next  X 

143:  Q-»Y 

144:  pit  5 , Y ,  1; cpl t  -2,-.3;lbl  Y; Y+500+Y ; if  Y*2Q00;gtc  146 

145:  gto  144 

146:  pit  3,-150 , 1; lbl  "Log  Intensity" ;plt  -1,500,1 

147:  csiz  1.7 ,1. 7, .71,90; lbl  "No.  or  Counts" ;csiz  1 .7 , 1. 7 , .71 ,0 

146:  for  J=1  tc  256 

149:  pit  a [ J ]  ,R[ J] *-2;  pen; next  J;pen;plt  0,0, -2 
150:  tor  J=1  to  256;plt  S[J]  ,P[J  ] ;  next  J 
151:  pit  6 . 2, 1169 , 1; lbl  “ .  Exper.  Curve" 

152:  pit  6.2,1122,l;lbl  "-  Cbi  Sq.  Prediction" 

153:  pit  8,1500,1 
154:  ret 
J55:  end 
*32056 
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APPENDIX  C 

CHI-SQUARE  GOODNESS  OF  FIT  TESTING  TO 

- THE  PI5TRIOTTTTO  HTP^THEST? - 

The  following  discussion  will  illustrate  the  rationale 
behind  our  selection  of  the  chi-square  test  to  examine  the 
assumed  log  normal  intensity  distributions  in  our  experiments. 
It  will  also  describe  the  methodology  utilized  in  the  actual 
determination  of  the  chi-square  value  that  was  calculated 
for  each  log  intensity  distribution.  Two  theorems  are 
necessary  in  order  to  begin.  [Ref.  7] . 

Theorem  1:  Assume  that  Y  has  distribution  function  F  , 
which  is  either  discrete  or  continuous,  and  that  Y1,  Y  , 

Yn  is  a  random  sample  of  Y.  Then  if  we  define 

Ix  *  (y:y  1  a  ),  I2  ■  Cy:a1  <y  <  a2)  ..., 

Ik-i  =  (y;ak-2  <y  -  ak-i)»  rk  *  Cy:ak-i  <  y)  » 

and  let  be  the  number  of  sample  values  that  fall  in  1^, 
i  3  l,2,....,k,  CXx ,  X2,  ...»  Xk)  is  the  multinomial  random 
variable  with  parameters  n,  p:,  p2,...,  p^,  where  n  is 
the  sample  size  (for  the  sample  of  Y)  and 

*  P[Y  falls  in  1^],  i  =  1,2,3, ...,k. 

Theorem  2:  Assume  that  Y  has  distribution  function  F 
which  satisfies  certain  regularity  conditions  and  which  is 
indexed  by  r  unknown  parameters,  0j,  02 ,  . 0  ,  and  that 
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Yj,  Y2,  . ..,  Yn  is  a  random  variable  of  Y.  Let  ©x,  0  ,..., 

A 

©r  be  the  maximum  likelihood  estimators  for  0j,  ©2,  . ..,  0  , 
respectively,  and  define  I : ,  I2,  ...»  1^  and  Xj,  X2,  ...» 
as  in  Theorem  1.  Then,  if  we  define 

?i  *  P [Y  falls  in  Ij,  i  =  1,2, ...,k 

where  01(  ©2,...,  ©r  replace  ©  l ,  ©2,...,  ©r  in  Fy,  the 
distribution  of  the  statistic 

V  =  k  (Xi  -  n  Pj)2 

E  n  P. 

i-1  1 

tends  to  that  of  a  chi-square  random  variable  with  k-r-1 
degrees  of  freedom  as  n  gets  large. 

Therefore,  when  it  is  desired  to  test  the  hypothesis 
that  Fy  has  a  specified  form,  the  actual  data  is  used  to 
estimate  the  parameters  and  then  the  chi-square  test  is 
initiated.  That  is,  the  observed  value  of  V  is  computed 
and  then  compared  with  the  100(1  -  a)  percentile  of  the 
chi-square  distribution  with  k-r-1  degrees  of  freedom.  If 
the  value  of  V  exceeds  this  latter  value,  then  the  hypo¬ 
thesis  that  Fy  has  the  specific  form  is  rejected. 

In  practice,  Larson  [Ref.  7]  states  that  there  is 
general  agreement  that  the  length  of  the  finite  intervals 
should  be  equal.  In  addition,  he  states  that  many  authors 
agree  that  the  chi-square  approximation  specified  in 
Theorem  2  is  quite  good  when  the  number  of  intervals  used  is 


as  large  as  possible,  consistent  with  the  constraint  that 
n  p^  ^  5  for  i«l,2,...,k. 

For  the  purposes  of  this  experiment,  it  was  desired  to 
test  the  hypothesis  that  the  log  intensity  distribution  for 
a  given  path  configuration  could  be  described  by  a  normal 
distribution.  Therefore,  the  hypothesis  was  that  the  log 
intensity  was  a  normal  random  variable. 

To  accomplish  this,  the  points  in  a  specific  signal 
distribution  (from  a  given  path  geometry  in  a  run)  were 
first  counted  then  used  to  compute  the  observed  values  of 
the  maximum  likelihood  estimators  for  the  normal  distribution. 
These  were  the  mean,  M,  and  the  variance,  a2,  and  they  were 
determined  as  follows: 

$  •  ~  ”  count  U) 

n 

i-  i 

and  a2  =  ^-  ^  (count  (z)  -  M)  2 

11  6 

1-  1 

where  n  ■  total  number  of  points  in 

a  distribution,  and 

count (z)  ■  intensity  value  of 
a  single  point. 

For  our  experiment,  and  subject  to  the  constraints 
discussed  previously  as  well  as  the  dynamic  range  of  the 
ADC  in  the  NIC-80,  256  intervals  of  equal  length  were  used 
with  n  *  16,384.  The  index  k  denotes  a  particular  interval 
or  bin.  Because  the  data  was  stored  in  a  histogram  of  counts 
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per  bin,  we  used  the  following  expressions,  which  are 
equivalent  to  those  listed  above,  to  compute  M  and  a2: 


M  =  i  2^6  (HKT)  Cfreq  00) 


n  z 
k=l 

1  256 

n  r 
k=l 


(freq  (K))  (bin  K  -  ft) 


where 


bin  k  =  mean  value  of  the  Kth  bin 
interval, 

freq(k)  =  number  of  counts  in  bin  K. 


Once  these  two  values  were  calculated,  they  were  used 
to  compute  the  observed  values  of  the  estimators  of  the 
probabilities  of  individual  log  intensity  values  falling 
into  the  256  bins  used  to  summarize  the  distribution.  To 
compute  these  values,  denoted  Px,  P2  ?256»  the 

probability  that  a  normal  random  variable  with  mean  and 

A 

standard  deviation  a  would  exist  in  each  of  the  256  intervals 
was  calculated. 

This  was  easily  accomplished  for  the  normal  distribution 
by  using  the  standard  normal  table  and  the  expressions 


Pk  .  IbxnJCh 


bin  kL  -  M 


where 


bin  kh  =  upper  value  for  bin  k, 
and 


bin  k^  =  lower  value  for  bin  k. 
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The  next  step  required  the  calculation  of  the  statistic  V, 
using  the  expression 

V  *  256  [freqOO  -  n(Pk)l2 
k-1  n(Pk) 

The  expression  for  degrees  of  freedom  from  Theorem  2  is 
degrees  of  freedom  3  k-r-1,  where  k  =  256,  the  number  of 

A  * 

bins,  and  r  =  2,  for  the  parameters  M  and  a  that  were 
estimated.  Therefore,  there  were  253  degrees  of  freedom 
for  this  process.  With  this  information,  the  final  step 
was  to  determine  the  value  for  the  desired  percentile  of 
the  chi-square  distribution  with  253  degrees  of  freedom  and 
compare  this  to  the  statistic  V.  We  chose  the  95th 
percentile,  with  a  value  of  287.88.  For  any  given  log 
intensity  distribution,  if  the  value  calculated  for  V  was 
less  than  287.88,  then  the  hypothesis  that  the  distribution 
was  normal  was  accepted.  If  V  was  greater  than  this  value, 
then  the  hypothesis  was  rejected. 


APPENDIX  D 


DATA  FROM  THE  FIRST  EXPERIMENT 


<R>(<1L_  >  ) 

S2 


Error (i 


— )  Runs 

N  Utilized 


Comments 


19 

All  sixteen  Two  splitters 

7 

All  (Decreasing 

ratios 

All  throughout 

the  run) 

6 

excl.  #3 

0 

All 

9 

All 

3 

All 

6 

excl.  #4 

1 

excl.  #6,12 

4 


NOV.  3  1982 

An, 1500, 2k, lk, 300k, 0 
N^ght,still,cir,outsd 
Wavelength*  6.33E-07  meters 


Scintillation  measurements  of  Sigma 
16384  samples  for  each  measurement 
Fixed  range*  100  meters 


I 


1 


s 


#  pts 

Sigma 

Chi  Square 

1:16385 

0.2321 

1.17E  03** 

1:16383 

0.1755 

8.92E  02** 

2:16385 

0.2109 

1.17E  03** 

2:16383 

0.1544 

1.13E  03** 

3:16385 

0.2159 

1.29E  03** 

3:16383 

0.1598 

1.34E  03** 

4:16385 

0.2162 

1.09E  03** 

4:16383 

0.1622 

9.98E  02** 

5:16385 

0.2162 

1.10E  03** 

5:16383 

0.1516 

1.17E  03** 

6:16385 

0.2154 

1.06E  03** 

6:16383 

0.1437 

1.56E  03** 

7:16385 

0.2043 

1.19E  03** 

7:16383 

0.1291 

1.37E  03** 

8:16385 

0.2206 

1.27E  03** 

8:16383 

0.1443 

1.27E  03** 

9:16385 

0.2017 

1.24E  03** 

9:16383 

0.1182 

1.36E  03** 

10:16385 

0.2038 

1  37E  03** 

10:16383 

0.1289 

1.59E  03** 

11:16385 

0.2076 

1.58E  03** 

11:16383 

0.1366 

1.49E  03** 

12:16383 

0.2070 

1.13E  03** 

12:16383 

0.1478 

1.38E  03** 

13:16383 

0.2015 

1.22E  03** 

13:16385 

0.1447 

1.27E  03** 

14:16383 

0.2057 

1.20E  03** 

14:16385 

0.1553 

1.07E  03** 

IS : 16383 

0.2094 

1.14E  03** 

15:16385 

0.1411 

1.29E  03** 

16:16382 

0.1957 

1.21E  03** 

16:16385 

0.1510 

1.70E  03** 

Figure  19.  Scintillation  Statistics 
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Figure  20.  Scintillation  Data  and  Theoretical  Curve 
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Scintillation  Data  and  Theoretical  Curve 


[C 


NOV.  2  1986 

Ap  > 1400 , 500 , lk, 300k, 10 
Night , cl r,slwd,outsd 
Wavelength*  6.33E-07  meters 


Scintillation  measurements  of  Sigma 
16384  samples  for  each  measurement 
Fixed  range*  100  meters 


#  pts 

Sigma 

Chi  Square 

1:16385 

0.1076 

7.99E 

02** 

1:16383 

0.1027 

7.67E 

02** 

2:16385 

0.1108 

9.7  7E 

02** 

2:16383 

0.1037 

1.12E 

03** 

3:16385 

0.1018 

1.23E 

03** 

3:16383 

0.1037 

7.62E 

02** 

4:16385 

0.1000 

7.80E 

02** 

4:16383 

0.1041 

7.19E 

02** 

5:16385 

0.1030 

5.98E 

02** 

5:16383 

0.0961 

8.17E 

02** 

6:16385 

0.1056 

5.10E 

02** 

6:16383 

0.0828 

1.55E 

03** 

7:16385 

0.1081 

4.89E 

02** 

7:16383 

0.0837 

1.24E 

03** 

8:16385 

0.1097 

5.14E 

02** 

8:16383 

0.1018 

8.33E 

02** 

9:16385 

0.1133 

6.19E 

02** 

9:16383 

0.1180 

6.2  3E 

02** 

10:16385 

0.1023 

8.66E 

02** 

10:16383 

0.1050 

7.54E 

02** 

11:16385 

0.0995 

5.96E 

02** 

11:16383 

0.0906 

1.25E 

03** 

12:16385 

0.0991 

7.87E 

02** 

12:16383 

0.0964 

1.35E 

03** 

13:16385 

0.1044 

8.08E 

02** 

13:16383 

0.1028 

9. 40E 

02** 

14:16385 

0.1101 

1.21E 

03** 

14:16383 

0.0978 

9.26E 

02** 

15:16385 

0.1077 

5.41E 

02** 

15:16383 

0.0924 

1.06E 

03** 

16:16385 

0.0993 

6. 25E 

02** 

16:16383 

0.0868 

1. 36E 

03** 

Figure  22.  Scintillation  Statistics 
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Theoretical  Curve 
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Figure  24.  Scintillation  Data  and  Theoretical  Curve 
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NOV.  3  1982 

Ap , 1500 ,G100 ,1k, 300k, 20 
Night , out sd,clr,wdy 
Wavelength*  6.33E-07  meters 


Scintillation  measurement  or  Sigma 
16384  samples  for  each  measurement 
Fixed  range*  100  meters 


#  pts 

Sigma 

Chi  Square 

1:16385 

0.0764 

4.18E 

02** 

1:16383 

0.0833 

6.27E 

02** 

2:16385 

0.0663 

5.36E 

02** 

2:16383 

0.0717 

1.03E 

03** 

3:16385 

0.0652 

4. 57E 

02** 

3:16383 

0.0710 

9.02E 

02** 

4:16385 

0.0642 

6.40E 

02** 

4:16383 

0.0689 

1.27E 

03** 

5:16385 

0.0614 

5.56E 

02** 

5:16383 

0.0603 

1.54E 

03** 

6:16385 

0.0751 

5.31E 

02** 

6:16383 

0.0719 

1.07E 

03** 

7:16385 

0.0636 

5.58E 

02** 

7:16383 

0.0683 

1.26E 

03** 

8:16385 

0.0609 

4.74E 

02** 

8:16383 

0.0674 

1.39E 

03** 

9:16385 

0.0824 

6.88E 

02** 

9:16383 

0.0781 

7.01E 

02** 

10:16385 

0.0800 

5.65E 

02** 

10:16383 

0.0787 

1.16E 

03** 

11:16385 

0.0829 

4. 8  3E 

02** 

11 : 16383 

0.0836 

9. 38E 

02** 

12:16385 

0^0979 

3.72E 

02** 

12:16383 

0.1015 

6. 53E 

02** 

13:16385 

0.1125 

2. 96E 

02** 

13:16383 

0.0936 

5.69E 

02** 

14:16385 

0.1116 

5.  32E 

02** 

14:16383 

0.1081 

5.90E 

02** 

15:16385 

0.1069 

5.48E 

02** 

15:16383 

0.0957 

6.14E 

02** 

16:16385 

0.1423 

5.97E 

02** 

16:16383 

0.1152 

5.  88E 

02** 

Figure  25.  Scintillation  Statistics 
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Figure  26.  Scintillation  Data  and  Theoretical  Curve 
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Figure  27.  Scintillation  Data  and  Theoretical  Curve 
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NOV.  3  1982 

AP, 1500, 500, lk, 10k, 30 
Night , still ,clr,outsd 
Wavelength-  6.33E-07  meters 


Scintillation  measurements  of  Sigma 
16384  samples  for  each  measurement 
Fixed  range*  100  meters 


#  pts 

Sigma 

Chi  Square 

1:16385 

0.1013 

6.24E 

02** 

1:16383 

0.0808 

8. 30E 

02** 

2:16385 

0.1022 

6.22E 

02** 

2:16383 

0.0799 

9.29E 

02** 

3:16385 

0.0997 

3.85E 

02** 

3:16383 

0.0700 

8.61E 

02** 

4:16385 

0.0948 

3. 86E 

02** 

4:16383 

0.0741 

7.99E 

02** 

5:16385 

0.0919 

6.22E 

02** 

5:16383 

0.0716 

8.67E 

02** 

5:16385 

0.1066 

4.26E 

02** 

6:16383 

0.0876 

8.92E 

02** 

7:16385 

0.1006 

3.74E 

02** 

7:16383 

0.0864 

1.18E 

03** 

8:16385 

0.1195 

5.S0E 

02** 

8:16383 

0.0918 

7. 36E 

02** 

9:16385 

0.1193 

5.62E 

02** 

9:16383 

0.0969 

1.64E 

03** 

10:16385 

0.1120 

4.34E 

02** 

10:16383 

0.0918 

9.05E 

02** 

11:16385 

0.1012 

4.69E 

02** 

11:16383 

0.0844 

1.21E 

03** 

12:16385 

0.0990 

4.78E 

02** 

12:16383 

0.0954 

1. 02E 

03** 

13:16385 

0.1163 

6.21E 

02** 

13:16383 

0.1018 

1.16E 

03** 

14:16385 

0.1104 

5.69E 

02** 

14:16383 

0.0877 

2.  71E 

03** 

15:16385 

0.0977 

1.18E 

03** 

15:16383 

0.0841 

2.61E 

03** 

16:16385 

0.0990 

6.18E 

02** 

16:16383 

0.0755 

5. 10E 

03** 

Figure  28.  Scintillation  Statistics 
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Figure  29.  Scintillation  Data  and  Theoretical  Curve 
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Figure  30.  Scintillation  Data  and  Theoretical  Curve 


I 


I! 

a 


NOV.  3  1982 

Ap , 1S00 , 200 , lk, 300k, 40 
Night, outsd,clr, st ill 
Wavelength*  6.33E-07  meters 

Scintillation  measurements  of  Sigma 
16384  samples  for  each  measurement 
Fixed  range*  100  meters 


#  pts 

Sigma 

Chi 

Square 

1:16383 

0.1219 

8.49E 

02** 

1:16385 

0.0806 

1.20E 

03** 

2:16383 

0.1260 

6.64E 

02** 

2:16385 

0.0956 

1.51E 

03** 

3:16383 

0.0962 

4. 88E 

02** 

3:16385 

0.0772 

1.57E 

03** 

4:16383 

0.0870 

6.01E 

02** 

4:16385 

0.0774 

1.29E 

03** 

5:16383 

0.1022 

4.40E 

02** 

5:16385 

0.0843 

1.21E 

03** 

6:16383 

0.1067 

5.64E 

02** 

6:16385 

0.0797 

1.67E 

03** 

7:16383 

0.0788 

5.18E 

02** 

7:16385 

0.0835 

9.98E 

02** 

8:16383 

0.0739 

7.22E 

02** 

8:16385 

0.0793 

1.00E 

03** 

9:16383 

0.0767 

8.63E 

02** 

9:16385 

0.0807 

2.  40E 

03** 

10:16383 

0.0800 

7.29E 

02** 

10:16385 

0.0756 

2.33E 

03** 

11:16383 

0.0677 

7.42E 

02** 

11:16385 

0.0651 

1.33E 

03** 

12:16383 

0.0673 

6.45E 

02** 

12:16385 

0.0687 

1.11E 

03** 

13:16383 

0.0623 

6.42E 

02** 

13:16385 

0.0635 

2.15E 

03** 

14:16383 

0.0657 

5.  59E 

02** 

14:16385 

0.0671 

1.98E 

03** 

15:16383 

0.0604 

1.  31E 

03** 

15:16385 

0.0679 

2. 86E 

03** 

16:16383 

0.0596 

6. 59E 

02** 

16:16385 

0.0667 

1.94E 

03** 

Figure  31.  Scintillation  Statistics 
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I'igure  32.  Scintillation  Data  and  Theoretical  Curve 


NOV.  3  1982 

Ap, 1500, SO, lk, 300k, 50 
Night, outsd,st ill, cir 
Wavelength*  6.33E-07  meters 


Scintillation  measurements  of  Sigma 
16384  samples  for  each  measurement 
Fixed  range*  100  meters 


ff  pts 

Sigma 

Chi  Square 

1:16384 

0.1048 

4.86E 

02** 

1:16384 

0.0673 

3.95E 

03** 

2:16385 

0.1145 

9.24E 

02** 

2:16383 

0.0898 

1.78E 

03** 

3:16385 

0.1110 

5.66E 

02** 

3:16383 

0.0956 

1.17E 

03** 

4:16385 

0.1098 

4.69E 

02** 

4:16383 

0.0944 

1. 35E 

03** 

5:16385 

0.1256 

5.01E 

02** 

5:16383 

0.1106 

6.41E 

02** 

6:16385 

0.1095 

6.6SE 

02** 

6:16383 

0.1077 

9.43E 

02** 

7:16385 

0.1066 

7.36E 

02** 

7:16383 

0.0949 

9.48E 

02** 

8:16385 

0.1045 

5.49E 

02** 

8:16383 

0.0691 

2.  56E 

03** 

9:16385 

0.1213 

7.46E 

02** 

9:16383 

0.1050 

9.07E 

02** 

10:16385 

0.0991 

5.59E 

02** 

10:16383 

0.0782 

1.47E 

03** 

11:16385 

0.1034 

9. 33E 

02** 

11:16383 

0.0775 

2.  59E 

03** 

12:16385 

0.1149 

5.72E 

02** 

12:16383 

0.0804 

1.63E 

03** 

13:16385 

0.1134 

6.  37E 

02*8 

13:16383 

0.0701 

2. 21E 

03** 

14:16385 

0.096 9 

8. 77E 

02** 

14 : 16383 

0.0602 

6. 83E 

03** 

15:16385 

0.0911 

8. 16E 

02** 

15:16383 

0.0621 

6. 61E 

03** 

16:16385 

0.0930 

7.44E 

02** 

16:16383 

0.0619 

7. 91E 

03** 

Figure  34.  Scintillation  Statistics 
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Figure  36.  Scintillation  Data  and  Theoretical  Curve 


NOV.  3  1982 

Ap, 1500, 200, lk, 300k, 60 
Night ,clr,outsd 
Wavelength*  6.33E-07  meters 


Scintillation  measurements  of  Sigma 
16384  samples  for  each  measurement 
Fixed  range*  100  meters 


If  pts 

Sigma 

Chi  Square 

1:16383 

0.1034 

4.96E  02** 

1:16385 

0.1007 

7.61E  02** 

2:16383 

0.1177 

5.42E  02** 

2:16385 

0.1146 

9.26E  02** 

3:16383 

0.1169 

4.47E  02** 

3:16385 

0.1174 

7.08E  02** 

4:16382 

0.1453 

4.47E  02** 

4:16385 

0.1434 

7.05E  02** 

5:16383 

0.1411 

4.67E  02** 

5:16385 

0.1369 

7.50E  02** 

6:16383 

0.1266 

4. HE  02** 

6:16385 

0.1287 

6.57E  02** 

7:16383 

0.1215 

4.64E  02** 

7:16385 

0.1140 

8.08E  02** 

8:16383 

0.1221 

6.55E  02** 

8:16385 

0.1188 

7.88E  02** 

9:16383 

0.1198 

5.09E  02** 

9:16385 

0.1119 

7.79E  02** 

10:16383 

0.1337 

4.95E  02** 

10:16385 

0.1291 

6.29E  02** 

11:16383 

0.1282 

5. 11E  02** 

11:16385 

0.1251 

7.84E  02** 

12:16383 

0.1296 

4.71E  02** 

12:16385 

0.1222 

1.03E  03** 

13:16383 

0.1276 

5.47E  02** 

13:16385 

0.1121 

8.39E  02** 

14:16383 

0.1027 

5.28E  02** 

14:16385 

0.0988 

9.65E  02** 

15:16383 

0.1085 

4.78E  02** 

15:16385 

0.0994 

1.00E  03** 

16:16383 

0.1236 

8.58E  02** 

16:16385 

0.1038 

1 . 03E  03** 

Figure  37.  Scintillation  Statistics 
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Figure  39.  Scintillation  Data  and  Theoretical  Curve 


NOV.  3  1982 

Ap, 1400, 100, lk, 300k, 70 
Night ,clr,outsd,slwd 
Wavelength*  6.33E-07  meters 


Scintillation  measurements  of  Sigma 
16384  samples  for  each  measurement 
Fixed  range*  100  meters 


#  pts 

Sigma 

Chi  Square 

1:16384 

0.1069 

6.42E  02** 

1:16384 

0.0729 

1.80E  03** 

2:16385 

0.1250 

6.46E  02** 

2:16383 

0.0816 

1.60E  03** 

3:16385 

0.1275 

6.70E  02** 

3:16383 

0.0857 

2.19E  03** 

4:16385 

0.1630 

6.86E  02** 

4:16383 

0.0921 

1.28E  03** 

5:16385 

0.1406 

5.83E  02** 

5:16383 

0.0895 

1.12E  03** 

6:16385 

0.1154 

7.4SE  02** 

6:16383 

0.0822 

2.40E  03** 

7:16385 

0.1105 

6.87E  02** 

7:16383 

0.0742 

5.07E  03** 

8:16385 

0.1087 

6.10E  02** 

8:16383 

0.0743 

4.18E  03** 

9:16385 

0.1201 

5.80E  02** 

9:16383 

0.0860 

2.14E  03** 

10:16385 

0.1592 

6.68E  02** 

10:16383 

0.1051 

9.76E  02** 

11:16385 

0.1457 

4.75E  02** 

11:16383 

0.1077 

8.63E  02** 

12:16385 

0.1396 

5.06E  02** 

12:16383 

0.1174 

7.79E  02** 

13:16385 

0.1435 

1.06E  03** 

13:16383 

0.0934 

1.45E  03** 

14:16385 

0.1238 

6.63E  02** 

14:16383 

0.0899 

1.48E  03** 

15:16385 

0.1200 

6.88E  02** 

15:16383 

0.0791 

2.64E  03** 

16:16385 

0.1147 

6.39E  02** 

16:16383 

0.0746 

3.99E  03** 

Figure  40.  Scintillation  Statistics 
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Scintillation  Data  and  Theoretical  Curve 
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Figure  42.  Sctintillation  Data  and  Theoretical  Curve 


NOV.  3  1982 

Ap, 1500, 200, lk, 300k, 80 
Night ,outsd, cl r,slwd 
Wavelength*  6.33E-97  meteres 


Scintillation  measurements  of  Sigma 
16384  samples  for  each  measurement 
Fixed  range*  100  meters 


#  pts 

Sigma 

Chi  Squar 

1:16384 

0.1045 

9.78E 

02** 

1:16384 

0.0849 

1.42E 

03** 

2:16383 

0.1175 

7.21E 

02** 

2:16385 

0.0952 

1.42E 

03** 

3:16383 

0.1011 

7.96E 

02** 

3:16385 

0.0914 

1.52E 

03** 

4:16384 

0.1170 

7.72E 

02** 

4:16384 

0.1049 

1.08E 

03** 

5:16385 

0.1206 

6.14E 

02** 

5:16383 

0.0942 

1.16E 

03** 

6:16385 

0.2736 

4.72E 

03** 

6:16383 

0.0920 

1.24E 

03** 

7:16385 

0.1003 

9.70E 

02** 

7:16383 

0.0828 

1.64E 

03** 

8:16385 

0.0931 

6.69E 

02** 

8:16383 

0.0762 

1.07E 

03** 

9:16385 

0.0735 

8.25E 

02** 

9:16383 

0.0651 

1.73E 

03** 

10:16385 

0.0745 

9. 36E 

02** 

10:16383 

0.0695 

1.  55E 

03** 

11:16385 

0.0830 

5. 55E 

02** 

11 : 16383 

0.0704 

1. 76E 

03** 

12:16385 

0.0977 

1.41E 

03** 

12:16383 

0.0688 

1.29E 

03** 

13:16383 

0.0818 

2.44E 

03** 

13:16385 

0.0838 

2. 64E 

03** 

14:16383 

0.0740 

4. 85E 

02** 

14:16385 

0.0731 

1.  37E 

03** 

15:16383 

0.0761 

5.42E 

02** 

15:16385 

0.0772 

1.27E 

03** 

16:16383 

0.0861 

4. 89E 

02** 

16 : 16  38S 

0.0794 

1.60E 

03** 

Figure  43.  Scintillation  Statistics 
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Scintillation  Data  and  Theoretical  Curve 
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Scintillation  Data  and  Theoretical  Curve 


NOV.  3  1982 

A,  1500, 200, 30k, 300k, 90 
Night ,outsd,clr ,slwd 
Wavelength*  6.33E-07  meters 


Scintillation  measurements  of  Sigma 
16384  samples  for  each  measurement 
Fixed  range®  100  meters 


1 


I 


#  pts 

Sigma 

Chi  Square 

1:16385 

Q. 1390 

5.99E 

02** 

1:16383 

0.1222 

6.37E 

02** 

2:16385 

0.1186 

8.47E 

02** 

2:16383 

0.1196 

6.72E 

02** 

3:16385 

0.0857 

8.55E 

02** 

3:16383 

0.0835 

1.02E 

03** 

4:16385 

0.0901 

7.97E 

02** 

4:16383 

0.0829 

1.22E 

03** 

5:16385 

0.0845 

6.44E 

02** 

5:16383 

0.0808 

1. 37E 

03** 

6:16385 

0.0863 

1.09E 

03** 

6:16383 

0.0860 

1 .  HE 

03** 

7:16385 

0.0900 

8.47E 

02** 

7:16383 

0.0861 

1.47E 

03** 

8:16385 

0.1117 

6.4  3E 

02** 

8:16383 

0.0808 

1.48E 

03** 

9:16385 

0.0938 

7.94E 

02** 

9:16383 

0.0800 

1.40E 

03** 

10:16385 

0.0961 

8. 37E 

02** 

10:16383 

0.0841 

1.59E 

03** 

11:16385 

0.1094 

8.  96E 

02** 

11:16383 

0.0987 

1.15E 

03** 

12:16385 

0.1046 

6.64E 

02** 

12:16383 

0.0966 

1.07E 

03** 

13:16385 

0.0983 

7.01E 

02** 

13:16383 

0.0897 

1. 38E 

03** 

14:16385 

0.0978 

8. 04E 

02** 

14:16383 

0.0814 

1.42E 

03** 

15:16385 

0.0959 

7.48E 

02** 

15:16383 

0.0869 

1.  5  3E 

03** 

16:16385 

0.1098 

6. 80E 

02** 

16:16383 

0.0984 

8. 24E 

02** 

Figure  46.  Scintillation 
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Scintillation  Data  and  Theoretical  Curve 


on  Data  and  Theoretical  Curve 


< 


NOV.  3  1982 

Ap, 1500, 100, lk, 300k, 100 
Night ,clr,outsd,slwd 
Wavelength®  6.33E-07  meters 


Scintillation  measurements  of  Sigma 
16384  samples  for  each  measurement 
Fixed  range®  100  meters 


#  pts 

Sigma 

Chi  Square 

1:16384 

0.0985 

9.13E 

02** 

1:16384 

0.1025 

2. 76E 

03** 

2:16385 

0.1019 

7. 58E 

02** 

2:16383 

0.1031 

2.04E 

03** 

3:16385 

0.1058 

9.34E 

02** 

3:16383 

0.1060 

1.94E 

03** 

4:16385 

0.1035 

9.79E 

02** 

4:16383 

0.1072 

2. 20E 

03** 

5:16385 

0.1004 

1.05E 

03** 

5:16383 

0.1076 

2.91E 

03** 

6:16385 

0.0966 

9.66E 

02** 

6:16383 

0.1043 

2.46E 

03** 

7:16385 

0.1103 

6.99E 

02** 

7:16383 

0.1158 

1.61E 

03** 

8:16385 

0.1083 

8.17E 

02** 

8:16383 

0.1126 

1. 81E 

03** 

9:16385 

0.1001 

1.05E 

03** 

9:16383 

0.1118 

2. 30E 

03** 

10:16385 

0.0949 

7.95E 

02** 

10:16383 

0.1115 

2. 35E 

03** 

11:16385 

0.0978 

8. 12E 

02** 

11:16383 

0.1136 

1.98E 

03** 

12:16385 

0.1198 

8.18E 

02** 

12:16383 

0.1276 

1.18E 

03** 

13:16385 

0.1181 

6.48E 

02** 

13:16383 

0.1263 

1. 36E 

03** 

14:16385 

0.1309 

6. 81E 

02** 

14:16383 

0.1397 

6. 57E 

02** 

15:16385 

0.1517 

4.01E 

02** 

15:16383 

0.1529 

6.75E 

02** 

16: 16385 

0.1537 

5.40E 

02** 

16:16383 

0.1441 

7. 99E 

02** 

Figure  49. 

Scintillation  Statistics 
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Scintillation  Data  and  Theoretical  Curve 
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DATA  FROM  THE  SECOND  EXPERIMENT 


D(cra) 

ap  ^ 

<R>  (<,s — 2r>) 
s 

Error (±  g) 

Runs 

Utilized 

5 

6.25 

.126 

All  sixteen 

10 

4.56 

.165 

excl.  #16 

20 

3.54 

.473 

All 

30 

1.90 

.098 

excl.  #8,  9, 
12,  13 

40 

4.32 

.096 

All 

50 

1.93 

.056 

excl.  #13,  14 
15,  16 

60 

5.16 

.167 

All 

70 

2.55 

.052 

All 

80 

5.87 

.144 

excl.  7 

90 

1.40 

.039 

All 

100 

4.98 

.088 

All 

110 

1.96 

.013 

All 

150 

6.01 

.121 

All 

126 


Comments 


Consistently 

decreasing 

ratios 

Four  pairs 
inverted  data 


Setup  mis¬ 
aligned  midrun 


i 


NOV.  13  1962 
Ap, 12O0G500, 100 ,5 
[  Nignt ,clr , roor 

wavelengtn*  6.33E-07  meters 


I 


1  4 


Scintillation  measurements  or  Sigma 
16384  samples  for  eacn  measurement 

40  meters 


fixed  range- 


#  pts 

S  igma 

Cni  Square 

1:16384 

0.0962 

3.24  E 

02** 

1:16384 

0.2507 

2.71E 

02 

2:16385 

0.0961 

4.48E 

02** 

2:16383 

0.2356 

4.40E 

02** 

3:16383 

0.0946 

2.67E 

02 

3:16378 

0.2497 

3.86E 

02** 

4:16383 

0.0983 

1.552 

02 

4:16385 

0.2469 

1.  7QE 

02 

5: 16382 

0.1021 

3.17E 

02** 

5:16379 

0.2675 

4.85E 

02** 

b:16383 

0.0885 

2.26E 

02 

6:16385 

0.2196 

2.60E 

02 

7:16383 

0.0931 

5.61E 

02** 

7:16385 

0.2359 

3.  756 

02** 

8:16383 

0.0988 

3.  72E 

02** 

8: 16385 

0.2432 

4.14E 

02** 

9: 16383 

0. 1119 

4.92E 

02** 

9:16385 

0.2661 

2.91E 

02** 

10:16382 

0.0981 

4.08E 

02** 

10:16382 

0.2470 

5.  76E 

02** 

11: 16383 

0.1107 

2.35E 

02 

11:16385 

0.2852 

2.25E 

02 

12:16383 

0. 1026 

3.  76E 

02** 

12:16385 

0.2558 

4.98E 

02** 

13:16384 

0.1122 

2 . 3o£ 

02 

13:  16363 

0. 2519 

4 . 05E 

02** 

14:16385 

0.0938 

3.  olE 

02** 

14:16383 

U.2284 

4.55E 

02** 

15:16384 

0.1063 

1.  90E 

02 

15 : 16381 

0.2704 

5. 35L 

02** 

16: 16382 

0.1206 

3.  10E 

02** 

16:16380 

0.2996 

3.65E 

02** 

Figure  52.  Scintillation  Statistics 


on  Data  and  Theoretical  Curve 


NOV.  13  19d2 
Ap, 1200G200 , 100 ,10 

Night , root  ,clr 

Wavelength-  6.33E-07  meters 


Scintillation  measurements  of  Sigma 
16384  samples  tor  eacn  measurement 
Fixed  range-  40  meters 


#  pts 

S  igma 

1:16310 

0.1401 

1:16090 

0.3073 

2:16333 

0.1348 

2:16160 

0.2802 

3:16328 

0.1191 

3:16165 

0.2694 

4:16331 

0.1199 

4:16189 

0.2475 

5: 16375 

0.0962 

5:16350 

0.2075 

6:16291 

0.1410 

6:16016 

0.2900 

7:16313 

0.1424 

7:16073 

0.2718 

8:16275 

0.1539 

8:15937 

0.2809 

9:16257 

0.1418 

9:15900 

0.2771 

10:16300 

0.1194 

10:16080 

0.2564 

11:16316 

0.1256 

11:16101 

0.2525 

12:16315 

0.1146 

12: 16110 

0. 26u2 

13:16361 

0.0802 

13:16294 

0.  1892 

14:16367 

0.0839 

14:16315 

0.1944 

15:16374 

0.0740 

15: 16344 

0.1716 

16:16368 

0 . 1 58  4 

16:16322 

0.1600 

Figure  55. 


Chi  square 

3.73E  02** 
1.S7E  02 
6.43E  02** 
2.67E  02 
4.62E  02** 
3.16E  02** 
8. 37  E  02** 
3.29B  02** 
3.26E  02** 
1.47E  02 
6.59E  02** 
3.44E  02** 
1.20E  03** 
2.34E  02 
1.82E  03** 
3.31E  02** 
1.02E  03** 
4.63E  02** 
4.16E  02** 
2.41E  02 
1.33E  03** 
4.58E  02** 
8. 99E  02** 
5.13E  02** 
4.63E  02** 
3.81E  02** 
7.33E  02** 
3.71E  02** 
3. 58E  02** 
3.81E  02** 
1.17E  04** 
3.aOE  02** 
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Figure  56.  Scintillation  Data  and  Theoretical  Curve 


NOV.  13  19b 2 

Ap, 1200, G50 0,100,20 

Nignt ,root ,clr 

wavelengtn=  b.33£-07  meters 


Scintillation  measurements  ot  Sigma 
16384  samples  tor  each  measurement 
fixed  range*  40  meters 


#  pts 

S  igma 

Cm  Squan 

1:16323 

0.1101 

4 . 88E 

02** 

1 :  lo200 

0.2694 

3.30E 

02** 

2:16353 

0.0790 

9.02E 

02** 

2:16319 

0.1996 

5.  29E 

02** 

3:16330 

0.0840 

6. 96E 

02** 

3:16263 

0 . 2 18  6 

1.63E 

02 

4:16310 

0.0814 

1.03E 

03** 

4:16238 

0.2095 

3. 74E 

u2** 

5:16153 

0.1545 

1.  70E 

03** 

5:15684 

0.3343 

4.28E 

02** 

6: 16218 

0.1199 

3. 55E 

03** 

6: 16021 

0.2426 

2.66E 

02 

7:16127 

0.1407 

4.97E 

03*  * 

7:15845 

0.2609 

3.65E 

02** 

8: 15894 

0.1904 

4.52E 

03** 

8:15210 

0.3360 

2.08E 

02 

9:15b91 

0. 1997 

4. 59E 

03** 

9: 14992 

0. 3312 

3. 13E 

02** 

10:15392 

0.2296 

4.14E 

03** 

10:14447 

0. 3601 

2.93E 

02** 

11:15118 

0.2508 

5.  35E 

03*  * 

11:14082 

0.364U 

3. 17E 

02** 

12:14630 

0.2961 

5.  98E 

03*  * 

12:13131 

0. 3  96  9 

5.  39L 

u2** 

13:13725 

0.3362 

6 . 3  4E 

03*  * 

13:11730 

0.4281 

1.05E 

03*  * 

14 : 15o94 

0.4231 

9 .  o  2b 

0  3*  * 

14:13359 

0.5351 

2 . 83E 

03*  * 

15  : 15  4  98 

u. 3319 

7.  39E 

03*  * 

15: 14U  30 

0.4365 

1 .  u9E 

o  3*  * 

16 : 15694 

0.2967 

1.  loE 

04*  * 

lo:i<*7d4 

0. 3523 

5.  2dc. 

02** 

Figure  58.  Scintillation  Statistics 
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Figure  S9.  Scintillation  Data  and  Theoretical  Curve 


I 


NOV.  13  1982 
Ap , 1200G1O0 , 3U 
Night,  root  ,clr 

wavelength*  6.33E-07  meters 


Scintillation  measurements  of  Sigma 
16384  samples  for  each  measurement 
Fixed  range*  40  meters 


1 


#  pts 

Sigma 

Chi  Squar 

1:16384 

0.1363 

5.25E 

02** 

1:16384 

0.2128 

3.47E 

02** 

2:16385 

0.1136 

5.  59E 

02** 

2:16383 

0 . 1 50  7 

7.69E 

02** 

3:16385 

0.0891 

1.08E 

03** 

3:16383 

0.1070 

2.69E 

03** 

4:16385 

0.1582 

1.65E 

03** 

4:16383 

0.2368 

3.  71E 

02** 

5:16385 

0.1122 

6.13E 

02** 

5:16383 

0.1612 

5.24E 

02** 

6:16385 

0.0937 

3.88E 

02** 

6:16383 

0.1067 

2.35E 

02 

7:16385 

0.1262 

9.15E 

02** 

7:16383 

0.1811 

7.92E 

02** 

8:16385 

0.0933 

4.47E 

02** 

8:16383 

0.0884 

3.20E 

02** 

9:16385 

0.1126 

2.57E 

03** 

9:16383 

0.0928 

1.18E 

03** 

10:16385 

0.1547 

1. 27E 

03** 

10:16383 

0.2256 

4.81E 

02** 

11:16385 

0.1424 

1.  91E 

03** 

11:16383 

0. 1712 

1.  02E 

03** 

12:16385 

0.1010 

6.20E 

02** 

12:16383 

0.0992 

1.05E 

02 

13:16385 

0.0980 

3.45E 

02** 

13:16383 

0.0962 

3. 01E 

u2** 

14:16385 

0.1145 

1.14E 

03** 

14:16383 

0.1483 

1.9  7E 

03*  * 

15:16385 

0.1461 

6. 05E 

02** 

15: 16383 

0.2085 

1.  42E 

02 

16:16385 

0.1631 

1.  UlE 

03** 

16 : 16383 

0.2432 

2. 43E 

02 

Figure  61.  Scintillation  Statistics 
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Scintillation  Data  and  Theoretical  Curve 


Scintillation  Data  and  Theoretical  Curve 


c 


NOV.  13  lid  2 
Ap, 120UG-T0U, 40 
Night ,  root , clr 

wavelengtn*  6.33t,-o7  n.eters 


Scintillation  measurements  or  Sigma 
16384  samples  ror  eacn  measurement 
Fixed  range*  40  n.eters 


#  pts 

Sigma 

Chi  Square 

1:16314 

0.1543 

2.60E 

02 

1:16132 

0.3395 

1.64E 

02 

2:16326 

0. 1470 

4. 93E 

02** 

2:16176 

0.3088 

1.64E 

02 

3:16353 

0.1227 

5.28E 

02** 

3:16279 

0.2618 

1.  99E 

02 

4:16331 

0. 1486 

3.86E 

02** 

4:16169 

0.3157 

2.33E 

02 

5:16323 

0.1528 

5.12E 

02** 

5: 16167 

0.3224 

2.58E 

02 

6:  16340 

0.1372 

3.40E 

02** 

6: 16209 

0.3083 

1.  78E 

02 

7:16362 

0.1219 

4.12E 

02** 

7:16292 

0.2756 

1.41E 

02 

8:16305 

0.1611 

2.39E 

02 

3:16088 

0.3302 

1.73E 

02 

9:16225 

0.1934 

4.24E 

02** 

9:15772 

0. 375S 

2.93E 

02** 

10:16323 

0. 1390 

8.  98  E 

02** 

10:16154 

0.2833 

5.  25E 

02** 

11:16249 

0.1357 

5. 54E 

02** 

11:15860 

0.3590 

2.12E 

02 

12:16275 

0.1705 

7.  85E 

02** 

12:15971 

0.3374 

1.  75E 

02 

13:16310 

0.1659 

6.  59E 

02** 

13:16080 

0.3384 

1.  82E 

02 

14:16258 

0. 1795 

5.9  4E 

02** 

14:15904 

0. 35e2 

2 . 64E 

02 

15:16305 

0.1751 

5 . 76E 

02** 

15: 1607b 

0.3410 

2 .  GiE 

02 

16:16371 

0.1137 

6 . 6  6E 

02** 

16: 16333 

0.2407 

4 . 25L 

02** 

Figure  64.  Scintillation  Statistics 
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l-'igure  66.  Scintillation  Data  and  Theoretical  Curve 


Ap,i2u4JlijS8,lo 

Night, clr, root 
Waveleng  tn* 


6.J3E-07  meters 


U 


Scintillation  measurements  of  Sigma 
16384  samples  for  each  measurement 
Fixed  range-  40  meters 


#  pts 

S  igma 

Chi  Square 

1:16278 

0.2057 

8. 56E 

02** 

1:16006 

0.2936 

3.49E 

02** 

2:16274 

0.2221 

7.15E 

02** 

2:15978 

0.3100 

1.61E 

02 

3:16224 

0.2467 

8. 58E 

02** 

3:15776 

0.3309 

2.34E 

02 

4:16293 

0.2219 

5.41E 

02** 

4:16032 

0. 3085 

2.01E 

02 

5:16282 

0.2217 

5.95E 

02** 

5: 16031 

0. 3243 

2.14E 

02 

6:16315 

0.1910 

7.92E 

02** 

6: 16135 

0.2893 

5.59E 

02** 

7:16337 

0.1988 

4.40E 

02** 

7:16194 

0.2907 

1.86E 

02 

8:16270 

0.2572 

1.43E 

03** 

8:15963 

0.3497 

4.39E 

02** 

9:16272 

0.2684 

1.45E 

03** 

9:15953 

0.3498 

2.23E 

02 

10:16214 

0.2790 

1.  30E 

03** 

10:15814 

0.3734 

2.54E 

02 

11:16082 

0.3507 

1. 86E 

03** 

11:15534 

0.4696 

3.  56E 

02** 

12:16274 

0.2541 

1.31E 

03** 

12:15989 

0.3367 

1.66E 

02 

13:  3851 

0.2811 

6. 27E 

02** 

13:12171 

1.7187 

6 . 88E 

U4  *  * 

Division  oy  zero  occurea  in  line  56 


Figure  67.  Scintillation  Statistics. 
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Figure  69.  Scintillation  Data  and  Theoretical  Curve 


NOV.  13  198  2 
Ap, 1200 ,G100  ,60 
Ntgnt, root ,clr 

waveleng tn»  6.33E-G7  meters 


Scintillation  measurements  or  Sigma 
16384  samples  tor  eacn  measurement 
Fixed  range-  40  meters 


#  pts 

S  igma 

Cni  ; 

Square 

1:16384 

0.1244 

3.  80E 

02** 

1:16381 

0.2827 

4.82E 

02** 

2:16385 

0.0984 

4.69E 

02** 

2:16383 

0.2254 

1.91E 

02 

3:16385 

0.0906 

5.  94E 

02** 

3:16383 

0.1967 

3.23E 

02** 

4:16385 

0.0626 

4.36E 

02** 

4:16383 

0. 1222 

2.85E 

02 

5:16385 

0.0712 

7.29E 

02** 

5:16383 

0. 1472 

3. 74E 

02** 

6:16385 

0.0941 

5. 10E 

02** 

6:16383 

0.2146 

6.82E 

02** 

7:16385 

0.1025 

5. 12E 

02** 

7:16383 

0.2289 

2.89E 

02** 

8:16385 

0.1091 

4.96E 

02** 

8:16383 

0.2486 

1.  75E 

02 

9:16385 

0.1083 

4.42E 

02** 

9:16382 

0.2442 

4.28E 

02** 

10:16385 

0.1203 

2.30E 

02 

10:16383 

0.2789 

2. 29E 

02 

11: 16385 

0.1071 

2.  52E 

02 

11:16383 

0.2556 

1.67E 

02 

12:16385 

0.0934 

2. 57E 

02 

12:16383 

0.2388 

2.61E 

02 

13:16385 

0.0947 

d .  6 6E 

02** 

13:16383 

0.2226 

2.93E 

02** 

14:16385 

0.0799 

8.  81 E 

02** 

14:16383 

0.1939 

7 .  72E 

02** 

15:16385 

0.0676 

b  .  20  E 

U2** 

15: 16383 

0. 1421 

1. 72E 

02 

16:16385 

0 . 0  7d9 

2.  55E 

02 

16:16383 

0. 1853 

2  .  U3E 

02 

Figure  70.  Scintillation  Statistics. 
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NOV.  13  1982 
Ap,12O0GlK,lUO,7u 
Night, root , clr 

wavelength®  6.33E-07  meters 


Scintillation  measurements  ot  Sigma 
16384  samples  tor  each  measurement 
Fixed  range*  40  meters 


#  pts 

S  igma 

Cni  Squar 

1:16364 

0.1687 

4.51E 

02** 

1: 16384 

0.2606 

2.48E 

02 

2:16385 

0.1517 

4 . 49E 

02** 

2:16383 

0.2360 

3.17E 

02** 

3:16385 

0.1266 

4.  73E 

02** 

3:16383 

0.1851 

5.  83E 

02** 

4:16385 

0.1169 

4 . 80E 

02** 

4:16383 

0.1543 

5.  59E 

02** 

5:16383 

0.1972 

6. 35E 

02** 

5:16382 

0.3320 

1.01E 

03** 

6:16380 

0.2272 

6.86E 

02** 

6: 16382 

0.3486 

4.10E 

02** 

7:16379 

0.2727 

9.05E 

02** 

7: 16373 

0.4092 

2.95E 

02** 

8:16383 

0.2343 

7.17E 

02** 

8:16383 

0.3525 

2.97E 

02** 

9:16385 

0.1706 

6.84E 

02** 

9:16383 

0.2398 

6.01E 

02** 

10:16385 

0.1601 

3.  53E 

02** 

10:16383 

0.2298 

2.9  3E 

02** 

11:16385 

0.1698 

7. 16E 

02** 

11: 16383 

0.2651 

4 . 65E 

02** 

12:16385 

0.1421 

3.  76E 

02** 

12:16383 

0.2187 

6. 85E 

02** 

13:16385 

0.1215 

3.  36E 

02** 

13:16383 

0.1669 

3.48E 

02** 

14:16385 

0.154  9 

3.  59E 

02** 

14:16383 

0.2210 

4.86E 

02** 

15:16385 

0.1519 

4  .  32E 

02** 

15:16383 

0.2376 

2 . 7t>E 

02 

16: 16385 

0 . 2U90 

6 .  3  IE 

02** 

16:16383 

0.  3179 

2  .  7 5E 

02 

Figure  73.  Scintillation  Statistics. 
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Scintillation  Data  and  Theoretical  Curve 


NOV.  13  198  2 
Ap, 1200 , GlUO  ,60 
Nignt.clr ,roof 

Wavelengtn*  6.33E-07  meters 


Scintillation  measurements  ot  Sigma 
16384  samples  tor  eacn  measurement 
Fixed  range*  40  meters 


#  pts 

Sigma 

Cni  Square 

1:16383 

0. 1202 

9.41E 

02** 

1:16381 

0.3218 

5.  36E 

02** 

2:16383 

0.1262 

4. 50E 

02** 

2:16385 

0.3341 

2.5  9E 

02 

3:16383 

0.1092 

6.34E 

02** 

3:16385 

0.2715 

3.60E 

02** 

4:16383 

0.0978 

7.57E 

02** 

4:16385 

0.2398 

6.38E 

02** 

5: 16383 

0.1434 

3.23E 

02** 

5:16385 

0.3518 

1.  96E 

02 

6:  16383 

0.1382 

3. 46E 

02** 

6:16385 

0.3391 

3.01E 

02** 

7:16383 

0.1273 

1.  44E 

03** 

7: 16385 

0.2316 

3.  86E 

02** 

8:16383 

0. 1140 

2.50E 

02 

8:16385 

0.2679 

2.69E 

02 

9: 16383 

0. 1080 

1. 55E 

02 

9: 16385 

0. 2626 

2.76E 

02 

10:16383 

0. 1054 

3.48E 

02** 

10:16385 

0.2528 

4.02E 

02** 

11:16383 

0.0918 

6 . 70E 

02** 

11:16385 

0.2176 

9. 67E 

02** 

12: 16382 

0. 1251 

1. 23E 

03** 

12 : 16380 

U . 30  4b 

1.  30 E 

03** 

13: 16385 

0.0871 

1.  97E 

02 

13: lo383 

0. 1986 

5.04E 

02** 

14:16385 

0.0997 

2. 82E 

02 

14:16383 

0.2231 

2. 23E 

02 

15:16385 

0.1145 

3  .  1 8  E 

0  4** 

15: 16383 

0.2650 

2  .  o  8  h 

02 

lb  : lb3  8  5 

0.1070 

2 . 4  9  E 

02 

io :  16  3  8  3 

0.2463 

J.b  7E 

02** 

Figure  7*.  Scintillation  Statistic 
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Scintillation  Data  and  Theoretical  Curve 


NOV.  13  19b  2 
Ap,1200G200,lK300K,90 
Nignt,  root  ,clr 

rtaveleng tn-  6.33E-07  meters 


Scintillation  measurements  or  Sigma 
16384  samples  for  eacn  measurement 
Fixed  range-  40  meters 


#  pts 

s  igma 

Chi  Squar 

1:16380 

0.1520 

4.6  2E 

02*  * 

1:16371 

0.1830 

6.6  9E 

02** 

2:16381 

0.1487 

5.  77E 

02** 

2:16377 

0.1698 

5.  34E 

02** 

3: 16382 

0.1725 

5.40E 

02** 

3:16375 

0.1922 

3.90E 

02** 

4:16378 

0 . 1 54  3 

3.17E 

02** 

4:16369 

0.1758 

2.46E 

02 

5: 16378 

0.1547 

4.35E 

02** 

5:16366 

0.1897 

6.  99E 

02** 

6:16365 

0.1976 

5.62E 

02** 

6: 16310 

0.2575 

1.46E 

02 

7:16378 

0.1725 

7.88E 

02** 

7: 16361 

0.2161 

3.83E 

02** 

8:16383 

0.1547 

6. 35E 

02** 

8:16383 

0.1779 

2.26E 

02 

9:16382 

0.1521 

3.81E 

02** 

9: 16373 

0.1705 

3.  71E 

02** 

10:16381 

0. 1437 

7.48E 

02** 

10:16370 

0.1691 

9.  89E 

02** 

11:16379 

0. 1813 

5.  55E 

02** 

11:16354 

0.2096 

2.78E 

02 

12:16383 

0.1133 

3.6  2E 

02** 

12:16385 

0.1289 

1.  56E 

0  3** 

13:16372 

0.1815 

1.09E 

03** 

13:16347 

0.2243 

fa.  38E 

02** 

14:16335 

0.2540 

1. 09E 

03** 

14:16216 

0.3055 

3 . 5  7  £ 

02** 

15:16367 

0.2289 

5.  8 7 E 

o  2*  * 

15: 16323 

0.2673 

1.45E 

U  2 

16:16373 

0.1919 

7. 13E 

0  2** 

16: lb344 

0.2322 

4 .  4  2E 

02 

Figure  79.  Scintillation  Statistics. 


Figure  80.  Scintillation  Data  and  Theoretical  Curve 
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Figure  81.  Scintillation  Data  and  Theoretical  Curve 


NOV.  13  19d2 
Ap, 1200G100 , 1UU 
Night , root ,clr 

Wavelength*  6.33E-Q7  meters 


Scintillation  measurements  ot  Sigma 
16384  samples  tor  eacn  measurement 
Fixed  range*  40  meters 


#  pts 

S  igma 

Chi  Squar 

1:16380 

0. 1349 

3.85E 

02** 

1:16371 

0.2994 

2.6d£ 

02 

2:16384 

0.1374 

3. HE 

02** 

2:16382 

0.2911 

2.35E 

02 

3: 16382 

0.  1484 

8.  26E 

02** 

3: 16379 

0.3248 

3.  76E 

02** 

4:16385 

0.1521 

2.64E 

02 

4:16383 

0. 3285 

1.95E 

02 

5:16384 

0.150  7 

5.95E 

02** 

5: 16382 

0.3289 

2.3  3E 

02 

6:16381 

0.1496 

5.  G9E 

02** 

6:16379 

0.3409 

2.35E 

02 

7:16382 

0.1194 

4. 9  IE 

02** 

7:16380 

0.2881 

2.44E 

02 

8:16383 

0. 1540 

7.92E 

02** 

8:16376 

0.3334 

1.  96E 

02 

9:16383 

0.1386 

4.  79E 

02** 

9: 16379 

0.2955 

2. 22E 

02 

10:16383 

0. 1475 

7.19E 

02** 

1U : 16379 

0.3283 

2.37E 

02 

11:16379 

0.1354 

d.  66E 

02** 

11:16372 

0.3178 

5.  98 E 

02** 

12:16375 

0.1743 

5.  39E 

02** 

12:16356 

0.3981 

2. 41E 

02 

13:16380 

0.1353 

3. 60E 

02** 

13:16375 

0.3218 

2.32k. 

02 

14:16380 

0.1575 

3. 21E 

02** 

14:16370 

0. 3484 

2 . 28 E 

02 

15:16369 

0.1860 

4. 53E 

02** 

15:16342 

0.4107 

2.44E 

02 

16:16383 

0. 1409 

3.  52E 

02** 

16: 16381 

0. 3038 

1.  95E 

02 

Figure  82.  Scintillation  Statistics. 
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Scintillation  Data  and  Theoretical  Curve 


log  iNrcNstrr 

Figure  84.  Scintillation  Data  and  Theoretical  Curve 
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NOV.  13  1982 

Ap, 1000 , G20U , 1K3U0K , 110 

Nignt,clr , root 

wavelengtn*  6.33E-07  meters 


Scintillation  measurements  of  Sigma 
16384  samples  for  eacn  measurement 
Fixed  range*  40  meters 


#  pts 

Sigma 

Cni  Squar< 

1:16379 

0.2205 

6. HE 

02** 

1:16360 

0.2979 

1.41E 

02 

2:16384 

0.1805 

3. 77E 

02** 

2:16376 

0.2543 

4.44E 

02** 

3:16360 

0.2719 

8. 33E 

02** 

3:16288 

0.3796 

7.15E 

02** 

4:16347 

0.2883 

4.81E 

02** 

4:16262 

0. 4213 

2.16E 

02 

5:16378 

0.2063 

4.3  3E 

02** 

5: 16362 

0.3164 

3.68E 

02** 

6:16363 

0.2332 

1.08E 

03** 

6: 16319 

0.3313 

7.67E 

02** 

7:16374 

0.2104 

4 . 556 

02** 

7:16356 

0.3154 

1.61E 

02 

8:16381 

0.2040 

6.81E 

02** 

8:16375 

0.2954 

4.38E 

02** 

9:16380 

0.1948 

4.12E 

02** 

9:16372 

0.2882 

3.03E 

02** 

10:16385 

0.1405 

8.98E 

02** 

10:16383 

0. 1725 

2.04E 

03** 

11:16384 

0.1533 

7.93E 

02** 

11:16382 

U. 1864 

1.  30E 

03*  * 

12:16383 

0.1610 

2.15E 

02 

12:16385 

0.2053 

1. 75E 

02 

13:16377 

0.2307 

b  .  4  7E 

02** 

13:16361 

0. 3240 

5.  99E 

02** 

14:16382 

0.2356 

3.  76E 

02** 

14:16371 

0. 3320 

2. 58E 

02 

15:16384 

0.2204 

5.  75E 

02** 

15:16379 

0.3132 

4  .  Q6E 

02** 

16:16376 

0.2471 

3.  58E 

02** 

16 : 1635b 

0. 3434 

1.67E 

u  /. 

Figure  85.  Scintillation  Statistics. 
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NOV.  13  1982 

Ap, 1200G500,100  ,150 

Nignt  ,cir ,  roor 

Wavelength*  6.33E-07  n.ecers 


Scintillation  measurements  of  Sigma 
16384  samples  for  eacn  measurement 
Fixed  ranges  40  meters 


#  pts 

S  igma 

Chi  Squar 

1:16384 

0.1387 

3.40E 

02** 

1:16384 

0.3606 

3.  73E 

02** 

2:16385 

0.1424 

2.12E 

02 

2:16383 

0.3727 

3. 53E 

02** 

3:16385 

0.1220 

1.48E 

02 

3:16383 

0.3060 

3. 33E 

02** 

4:16385 

0.0941 

1.39E 

02 

4:16383 

0.2423 

4. HE 

02** 

5:16385 

0.1020 

1.  97E 

02 

5: 16383 

0.2403 

5.47E 

02** 

6:16385 

0.1165 

3.59E 

02** 

6:16383 

0.2799 

4.82E 

02** 

7:16385 

0.1241 

4.67E 

02** 

7:16383 

0.2895 

6.26E 

02** 

8:16384 

0.1322 

6.34E 

02** 

3:16382 

0.3146 

5.49E 

02** 

9:16383 

0.1322 

1. 78E 

02 

9: 16385 

0.3226 

4.91E 

02** 

10: 16383 

0. 1277 

2.45E 

02 

10:16385 

0.3108 

2.20E 

02 

11: 16383 

0.1090 

1.  96E 

02 

11:16385 

0.2531 

3. 20E 

02** 

12:16383 

0. 1237 

1.  91E 

02 

12:16385 

0.3059 

3. 03E 

02** 

13:16383 

0.0993 

2.21E 

02 

13: 16385 

0.2504 

1.  87E 

02 

14:16383 

0.1011 

2. 80E 

02 

14:16385 

0.2466 

3. 97E 

u  2** 

15: 16383 

0.0972 

3.  40E 

U2** 

15:16385 

0.2360 

2.64E 

02 

16 : 16383 

0. 0*91 

7 .  *  4E 

02** 

16:16385 

0.2382 

4  .  IdE 

02** 

Figure  38.  Scintillat 


Figure  89.  Scintillation  Data  and  Theoretical  Curve 
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Figure  90.  Sc int i 1 lat ion  Data  and  Theoretical  Curve 


APPENDIX  F 


DATA  FROM  THE  THIRD  EXPERIMENT 


<R>  y>) 

aS 

Error (+  — ) 

N 

Runs  Utilized 

2.01 

.035 

All  four 

2.46 

.073 

All 

1.65 

.018 

All 

2.24 

.091 

All 

1.44 

.050 

All 
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NO  V  .  21  l9b2 

Ap  ,  12u0  ,(j1k,  100,5 
Nignt  ,ovc , rool , stl 
wavelength*  6.33E-07  meters 


Scintillation  measurements  or  Siyir.a 
163b4  samples  tor  eacn  measurement 
Eixea  range*  4u  meters 


*  pts 

Sigma 

Chi  j 

Squar 

1:15842 

0.2071 

3.69E 

02** 

1:15848 

0.2906 

2.90  £ 

02** 

2:15799 

0. 2041 

3.  86E 

02** 

2: 15810 

U .  2979 

5.4  5E 

02** 

3:15816 

0.2251 

1.84E 

02 

3:15830 

0.3155 

1.  o9E 

02 

4:15831 

0.2105 

2. 90E 

02** 

4:15831 

0.2961 

1. 46E 

02 

Figure  91.  Scintillation  Statistics 


Figure  92.  Scintillation  Data  and  Theoretical  Curve 
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Figure  93.  Scintillation  Data  and  Theoretical  Curve 


NOv.  21  1962 
Ap,120G,Glk, luO,10 
Mgnt  ,cvc ,  root,  stl 
wavelengtn*  6.33E-o7  meters 


Scintillation  measurements  or  Sigma 
16364  samples  tor  eacn  measurement 
Eixea  ranges  40  meters 


#  pts 

Sigma 

Cm  Square 

1:15609 

0. 2038 

3.  71E 

02** 

1:15604 

0.3226 

2.61E 

02 

2:15816 

0.1822 

1.  46E 

02 

2:15817 

0.2865 

1.  78E 

02 

3:15844 

0.2249 

4.43E 

02** 

3:15646 

0.3359 

1.6  9E 

02 

4:15806 

0.1927 

3.  08E 

02** 

4:15808 

0.3128 

2.  77E 

02 

Figure  94.  Scintillation  Statistics 
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NCV.  21  19*2 
Ap,12U0,Glk, 100,20 
N ignt  ,ovc, root ,stl 
wavelengtn*  6.33E-07  meters 


Scintillation  measurements  ot  Sigma 
16364  samples  tor  eacn  measurement 
Eixea  ranges  40  meters 


f  pts 

Sigma 

Chi  squar 

1:1565b 

0.1767 

3.  77E 

02** 

1:15857 

0.2250 

5.65E 

02** 

2:15646 

0.1794 

2.43E 

02 

2:15849 

0.2285 

3.  08E 

02** 

3:15854 

0.1630 

5.23E 

02** 

3:15856 

0.2093 

4.03E 

02** 

4:15843 

0.1812 

2.89E 

02** 

4:15844 

0.2367 

4.24E 

02** 

Figure  97.  Scintillation  Statistics 
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Figure  99.  Scintillation  Data  and  Theoretical  Curve 


NOV.  21  19fa2 
Ap,12u0,Glk,100,30 
N  ignt ,ovc, root,  stl 
Wavelength*  6.33E-G7  meters 


Scintillation  measurements  or  Sigma 
16384  samples  tor  each  measurement 
P.ixea  range*  40  meters 


I  pts 

1:15834 

1:15835 

2:15850 

2:15856 

3:15847 

3:15846 

4:15841 

4:15841 


Sigma 

0.1801 
0.2612 
0. 20  50 
0.3195 
0.2064 
U. 3202 
0.2070 
0.2943 


Chi  Square 

2.35E  02 
2.96E  02** 
1.97E  02 
4.55E  02** 
2.02E  02 
2.27E  02 
1.83E  02 
2.56E  02 


Figure  100.  Scintillation  Statistics 
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Figure  102.  Scintillation  Data  and  Theoretical  Curve 


NOV.  21  19b2 
Ap,12UU ,Glk, 100, 40 


Mignt  ,ovc ,  root 

,  stl 

waveleng  tn= 

6 .33E-07 

meters 

Scintillation 

measurements  or 

Sign,, 

16384  samples 

tor  eacn 

measurement 

Fixea  ranges 

40  meters 

f  pts 

Sigma 

Chi  Sguar 

1:15851 

0.2509 

3.62E 

02** 

1:15850 

0.2859 

1.  97E 

02 

2:15o55 

0.2383 

4. 22E 

02** 

2:15854 

0.2826 

1.  53E 

02 

3:15851 

0.2208 

2.4  2E 

02 

3:15851 

0.2773 

3.13E 

02** 

4:15851 

0.2244 

8.09E 

02** 

4:15853 

0.2703 

6 . 08E 

02** 

Figure  103.  Scintillation  Statistics 
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Figure  104.  Scintillation  Data  and  Theoretical  Curve 


Figure  105.  Scintillation  Data  and  Theoretical  Curve 
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